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abstract 

The quality of semiconductor crystals grown by the vertical Bridgman 
technique is strongly influenced by the axial and radial variations of 
temperature within the charge. The present work examines the relationship 
between the thermal parameters of the vertical Bridgman system and the 
thermal behavior of the charge. 

Thermal models are developer which are capable of producing results 
expressable in analytical form and which can be used without recourse to 
extensive computer work for the preliminary thermal design of vertical 
Bridgman crystal growth syste-s. These models include the effects of 
thermal coupling between the larnace and the charge, charge translation 
rate, charge diameter, thickness and thermal conductivity of the confining 
crucible, thermal conductivity change and liberation of latent heat at the 
growth interface, and infinite charge length. The hot and cold zone 
regions, considered to be at spatially uniform temperatures, are separated 
by a "gradient control" region which provides added thermal design 
flexibility for controlling the temperature variations near the growth 
interface. 

One-dimensional thermal models show that the thermal coupling between the 
furnace and the charge and the change in thermal conductivity at the growth 
interface are the primary factors influencing the axial temperature 
gradients near the interface. The effect of the crucible is analytically 
described in a way which does not compromise the simplicity of the one- 
dimensional models; its effect is to reduce the thermal coupling between 
the charge and the fu’.nace. Criteria for the neglect of charge translation 
rate, charge laigth and liberation of latent heat are developed. 

Two-dimensi(xial models show that the interface shape is primarily affected 
by thermal conditions near the growth interface. Heat transfer to or from 
the charge near the interface is shown to provide a means for altering the 
interface shape. The combined effects of a crucible and the liberation of 
latent heat and/or conductivity change at the growth interface produce 
radial temperature variations at the growth interface. For the growth of 
semiconductors, this effect leads to an adverse interface curvature and 
remains an unresolved problem inherent to the Bridgman technique. 
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cross sectional area; without subscript denotes the charge cross 
sectional area 
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eq. [6.1]) 
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thermal diffusivity; without subscript denotes the thermal 
diffusivity of the charge 

Dqj-/D 

surface emissivity 
dimensionless aucial coordinate, Z/D 
axial location of the cold end of the charge 
axial location of the hot end of the charge 
axial location of the growth interface 
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The basic material property requirements of semiconductors used for 
advanced electronic devices are; 

(1) absence of grain and twin boundaries (single crystal) 

(2) low crystal defect density 

(3) uniform distribution of the dopant 

(4) minimized contamination by electrically active elements 

Numerous techniques exist for the productiai of semiconductors in single 
crystal form (c.f., [1,2]). The choice of a technique for a particular 
aKJli'ation is dictated b^- the inherent physical and chemical properties of 
the sei’iiconductor and the desired resultant electronic properties [3]. 
Increasingly stringent requirements for the electronic properties of 
semiconductor crystals necessitate the modification and refinement of 
conventional growth techniques. 

Growth from the melt describes several techniques used for the pro- 
duction of bulk crystals by the controlled transformation of a liquid into 
a solid of high degree of crystalline perfection. The majority of bulk, 
single crystal semiconductors is currently grown from the melt. The 
commonly used melt growth techniques are depicted in Fig. 1.1. In each 
case, growth of a single crystal is initiated by partial melting of a seed 
crystal; continued growth of the crystal occurs by the relative motion of 
the charge (crystal and melt) through the thermal boundary conditions 
imposed by the apparatus. 

Each of the melt-growth techniques has advantages and shortcomings. 
The Czcchralski technique (Fig. 1.1a), for example, has relatively high 
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rates of production and can produce crystals with a low defect density. 
Axial and radial segregation of the dopant is typically significant, 
however, due to the effects of time-varying, natural convective fluid flow 
in the melt. Also, precaution must be taken to prevent depletion of 
volatile elemeits from the melt. 

Crystals grown by the float-zone technique (Fig. 1.1b) are relatively 
free from contamination since the melt is uncontained; radial segregatic^ 
of the dopant is appreciable, however, due to the strong two-dimensional 
thermal field near the growth interface and the resulting natural 
convective fluid flow. 

Use of a sealed crucible in the vertical Bridgman technique 
(Fig. 1.1c) provides a simple means of minimizing the depletion of volatile 
elements from the melt. (The crucible used in the horizontal Bridgman 
technique is typically open.) Crystals grown by the Bridgman techique, 
however, have high defect densities generally attributed to mechanical 
stresses associated with the confining crucible. The boundary between the 
growth interface and the crucible also provides sites for heterogeneous 
nucleation of crystal grains which can prevent single crystal growth. 
Further, although the vertical Bridgman configuration is relatively stable 
for natural convective fluid flow (the melt is heated from above), radial 
segregation of the dopant in Bridgman grown crystals has not been 
demonstrated to be consistently superior to that in crystals produced by 
other melt-growth techniques. 

It is evident that heat transfer is involved in the phase transforma- 
tion at the crystal/melt (growth) interface and it is now understood that 
thermal phenomena greatly affect the quality of melt-grown crystals. Basic 
thermal considerations have not, however, been adequately applied in the 
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Fig. 1.1b ; Sketches of the ccrmonlv used rnelt graArth 
techniques : float-cone . 


- 23 - 


OF POOR QUALITY 



Fig. 1.1c; Sketches of the coirmonl'^ used nnelt gra-rth cechnicnaes 
vertical Bridgran. 
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analysis or design of melt-growth apparatus. There is, therefore, a strong 
motivation for a reassessment of the inherent capabilities of the various 
melt-growth techniques from the point of view of improved thermal design. 

Compared to other melt-growth techniques, the thermal configuration of 
the vertical Bridgman technique offers several features that appear 
advantageous for semiconductor crystal growth: 

(1) The thermal boundary conditions associated witli the furnace 
are relatively simple and axisymmetric. 

(2) The furnace configuration is readily adaptable to thermal 
design changes deemed necessary to accommodate the growth of 
different semiconductor materials with varying thermal 
properties. 

(3) Natural thermal convection in the melt is reduced since the 
melt is heated from above and is completely enclosed by the 
crucible. 

For these reasons a combined analytical and experimental program is being 
conducted to investigate the potential for improvement in the thermal 
design of the vertical Bridgman crystal growth technique. The present work 
is the analytical portion of this program. The methodology of the analysis 
if to develop simplified thermal models which are capable of showing the 
governinq parametric behavior of the growth system. Particular emphasis is 
given to obtaining approximate analytical expressions which easily 
demorstrate the interaction among the thermal parameters. 

Thermal requirements for satisfactory melt-growth are shown in section 
2.1 to depend on both axial and radial temperature variations within the 
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charge. One-dlmensicnal (i.e., axial tenifc?erature variations only) and two- 
dimensional thermal models are formulated in Chapters 4 and 6 respectively. 
The temperature field within the charge obtained from the models is th«i 
used to evaluate the effects of the thermal parameters from the viewpoint 
of satisfying the thermal requirements (Chapters 5 and 7). 
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2^ SOnPF. OP TOE TOERMAL ANALYSIS 

A combined analytical and experimental proqram for the improvement of 
semiconductor crystal growth by the vertical Bridgman technique is oeing 
undertaken in the Department of Materials Science and Engineering. The 
analysis presented in this thesis comprises the analytical portion of this 
program. 

The objective of the thermal analysis is to identify and provide an 
understemding of the thermal characteristics of semiconductor crystal 
growth by the vertical Bridgman technique; in this way, problems likely to 
be encountered, and the means for their possible reslolution, are expected 
to be identified. To perform this task requires, of course, an adequate 
knowledge or the effects of the principal thermal parameters. Models are 
formulated to ascertain these effects with special emphasis directed 
towards obtaining results which cotn be presented in analytical form and can 
th< refore be used without recourse to extensive computer work. Design 
alternatives ‘Jv’ggested by the tnermal analysis are to be evaluated in the 
experimental program. 

Section 2.1 points out that the benavior of the axial and radial 
temperature gradients near the growth interface is of primary importance in 
determining the quality of crystals grown by the Bridvman technique. The 
results of the thermal analysis are therefore discussed in this context. 
Chapters 4 and 5 deal ’'ith the axial temperature distributiai of ti'.e charge 
and Chapters 6 and 7 treat radial temperature variations within the charge. 
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2J. THERMAL REQUII^JimS FOR BRIDGMAN CRYSTAL GRCWTH 

The electronic properties of a crystal grown by a melt-growth 
technique depend strongly upon the temperature distribution near the growth 
interface during solidification. The thermal requirements leading to the 
growth of high quality semiconductor crystals by the vertical Bridgman 
technique are: (1) that the axial temperature gradient at the growth 
interface is greater than a critical value and (2) that the radial 
temperature gradients near the growth interface produce an interface shape 
which is slightly concave toward the solid (Fig. 2J.). 

2.1.1 Axial Temperature Gradient At The Growth Interface 

During solidification, a concentration boundary layer of the dopant is 
generated in the melt in advance of the growth interface. Within this 
boundary layer the solidification temperature varies in accordance with the 
phase diagram. The phenomenon termed constitutional supercooling occurs 
when the sensible temperature in the melt is less than the local 
solidification temperature (4]. Under such conditions the morphology of the 
interface becomes unstable and it is difficult to maintain single crystal 
growth 15,6,7]. If the concentration boundary layer in the melt is 
diffusion dominated (i.e., neglecting natural convective fluid flow), the 
axial temperature gradient in the melt at the interface required to prevent 
constitutional supercooling is [81: 



where: R = growth rate 

mL = slope of the liquidus of the phase diagram 
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where: Cq ■ concentration of dopant in the nvelt 

Kq ■ equilibrium partition ratio of the phase diagram 
d ■ diffusion coefficient of the dopant in the i!ielt 

Equation [2.1] indicates that the axial gradient criteria is, to a certain 
extent, controlled by the crystal grower. By decreasing the. lowering rate, 
V, the growth rate, R, likewise decreases [9,10,111 and the axial gradient 
criteria becomes less restrictive. 

llie axial gradient criteria is less restrictive in the presence of a 
convective compaient to the mass transfer of dopant, near the interface [4]. 
The axial gradient requirement given by eq. [2.11 represents, therefore, an 
upper limit. 

2.1,2 Shape of the Growth Interface 

Heterogeneous nucleation can occur at the boundary betwe«i the growth 
interface and the crucible. When allowed to grow, crystal grains initiated 
by such nucleation destroy single crystal growth. These unwanted crystal 
grains do not propagate into the single crystal, however, if the interface 
is concave toward the crystal as shown in Fig. 2.1. It is desired, 
therefore, that radial temperature variations within the charge near the 
growth interface produce an interface which is slightly concave toward the 
solid. 

The shape of the interface is also a factor concerning the radial 
segregation of the dopant in the solidified charge. For the special case 
of a diffusion dominated concentration boundary layer at the growth 
interface, analysis shows [12,131 that the radial mass transfer of the 
dopant decreases with interface curvature. Thus, to achieve a desired 
level of radial segregation under such conditions, it is only required that 
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th« interface curvature be reduced to an appropriate level. For dilutely 
doped semiconductors^ even in the presence of radial segregation, the 
interface is very nearly an isotherm. The interface shape criteria 
therefore translates into a thermal criteria — the interface isotherm 
should be nearly flat. That is, heat transfer near the interface should be 
marly one-dimensional. 

Cases of practical significance occur in which the assumptions leading 
to the flat interface criteria discussed above are not valid. For example, 
the theoretical studies of Chang [141 and Chang and Brown [151 show that 
natural convective fluid flows may have a significant effect on the radial 
segregation of the dopant while having only a minor effect on the 
temperature distribution in the charge for melts which have small Prandtl 
numbers. Under these conditions a flat growth interface is not directly 
correlated with satisfactory radial segregation of the dopant. Also, in 
systems which exhibit a large separation of the liquidus and solidus (e.g., 
HgCdTe) , the interface shape is coupled to the radial segregation of the 
dopant as well as the radial variation of the temperature in the vicinity 
of the interface. Such cases as these require a detailed analysis of the 
mass, momentum and heat transfer within the charge in order to predict the 
radial segregation. 

Dopant mass transfer is not within the scope of the present thermal 
analysis. The behavior of isotherm shapes to be described in subsequent 
chapters may therefore be correlated to radial segregation behavior only 
for the special case of dilutely doped, diffusion dominated systems. When 
these ideal conditions are not approximated, results from the present 
analysis are nevertheless expected to describe the thermal behavior '>f the 
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charge and may thus be used in conjunction with analyses which consider 
mass transfer of the dopant (e.g., [14,15]) to aid in understanding the 
manner in which radial segregation is affected by the thermal parameters of 
the Bridgman growth system. 


2.2 THE ■] 




IftL MODEL 


The two-dimensional temperature field in the solidifying charge is 
described by partial differential equations subject to complex boundary 
conditions. The solution of these equations may be effected by standard 
numerical techniques (e.g., finite difference or finite element methods). 
Numerical solution methods preclude, however, the description of results in 
simple analytical form. The present work emphasizes the development of 
approximate models which lead to analytical results. The increased 
accuracy of more exact numerical formulations is not considered warranted 
at the present time for the following reasons: 


(1) Thermal property values (e.g., thermal conductivity of the 
melt, emissivities of the furnace wall and crucible) are 
often not well known and introduce uncertainty into the 
results of any thermal analysis. 

(2) The development of approximate thermal models itself leads to 
a better understanding cf the thermal behavior of the 
Bridgman system. 

(3) Approximate results presented in analytical form and 
indicating the primary thermal behavior of the Bridgman 
system are more conveniently used to evaluate the effects of 
various furnace design alternatives. 
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The present work models the vertical Bridgman configuration being 
constructed for the experimental portion of the program. The experimental 
apparatus, described in section 2.2.1, contains the features basic to all 
vertical Bridgman growth systems; therefore, results obtained in this work 
can be used in a qualitative manner for any vertical Bridgman furnace 
configuratiOT. It is necessary, of course, that simplifying approximations 
be made in the development of any thermal model; the thermal effects 
incorporated in and excluded from the present models are outlined in 
Section 2.2.2. 


2.2.1 The Experimental Vertical Bridoman Crystal Growth System 

Vertical Bridgman crystal growth denotes the solidification of a 
molten charge, contained in a chemically inert crucible, by its descent 
from a hot to a cold environment. In conventional Bridgman growth 
(Fig. 1.1c) the hot environment is a tube furnace of varying design and the 
cold environment is the surroundings (e.g. the laboratory). The present 
work formulates a thermal model of a system, currently under construction, 
which incorporates several modifications of a conventional Bridgman system. 
There are two objectives for these modifications: 

(1) To permit flexibility for controlling the thermal environment 
exposed to the charge. In this way, through the adjustment 
of design and operation parameters, the system is expected to 
provide for a wide range of axial temperature gradients, and 
radial temperature gradients which lead to interface shapes 
ranging from convex to concave. 

(2) To permit the formulation of more accurate thermal models of 
the system, particularly in describing the thermal boundary 
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conditions exposed to the charge. In this way, the thermal 
analysis is expected to provide more applicable results. 

The experimental Bridgman system is depicted in Fig. 2.2. The cold 
aivironmoit is a separate furnace whose temperature can be independeitly 
controlled. Stockbarger [16] first introduced a separate furnace for the 
cold environment which was separated from the hot furnace by a simple 
radiatioi shield as a means for increasing the axial temperature gradient 
at the growth interface. In the present system, heat pipes are used for 
the hot and cold furnaces. The temperature boundary condition is therefore 
uniform in both the axial and circumferential directions. The gap between 
the heat pipes and the charge can be made smaller in order to simplify the 
description of radiant heat transfer and also, if desired, to obtain 
significant heat transfer by thermal conduction. Between the heat pipes is 
a "gradiait control" region in which the growth interface is to be located. 
Chang and Wilcox [17] suggested the use of an adiabatic gradient control 
region as a means of reducing the curvature of isotherms near the growth 
interface. The present work considers the thermal design of the gradient 
control region as a design variable. It is anticipated that judicious 
design of heat transfer between the gradient control region and the charge 
will assist in obtaining the requisite control over the interface shape. 

2.2.2 Effects Incorporated In and Excluded Frcm the Thermal Models 

The following thermal parameters, considered to be of primary 
importance in determining the thermal characteristics of Bridgman crystal 
growth, are included in the thermal models of Chapters 4 and 6: 
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Fig, 2.2 ; The experirental \’ertical Bridcanan cr'-stal 
gra\Tth s^'stem. 
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(1) Thernal coupling between the charge and the 
heat pipe furnace. 

Thermal communication between the charge and the 
furnace is due to both radiation and conduction heat 
transfer. The radiative contribution strongly depends 
on the values of surface emissivities which are 
oftentimes not accurately known. The conductive 
contribution is determined by the thickness of the gap 
between the charge and furnace and the conductivity of 
the gas in the gap. ThuSr evacuating the gas from the 
gap can eliminate the conduction heat transfer. 
Appendix C discusses the quantitative evaluation of 
typical heat transfer coefficients associated with both 
modes of heat transfer across the gap. 

(2) Thermal design of the gradient zone region. 

The gradient zone portion of the furnace is modelled as 
an annular ring between the hot and cold heat pipes. 
Changes in the thermal conductivity and thickness of 
the annular ring are expected to alter the distribution 
of heat transfer to the charge and, as a consequence, 
the shape of the isotherms within the gradient zone. 

(3) Thermal conductivity of the charge and the 
ratio, Rk, between the melt and crystal 
conductivities . 

For semiconductors, the melt conductivity is greater 
than that in the solid, i.e., Rr > 1. 

(4) Generation of latent heat at the growth 
interface. 
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(5) Thickness and thermal crnductivlty of the 
crucible. 

(6) Lowering rate of the charge. 

Lowering rates for typical semiconductor crystal growth 
range from 0,1 m/s to lO^m/s. 

(7) Length of the charge and its position in the 
furnace. 

(8) Diameter of the charge. 

Several effects have not been included in the thermal models. The 
excluded effects are of relatively small importance for typical crystal 
growth situations or do not contritxjte to the primary thermal behavior of 
the Bridgman growth system; it is therefore considered that the increased 
complexity of thermal models that would be required for inclusion of these 
effects is not warranted. The excluded effects are: 

(1) Natural convection in the melt. 

The temperature field in the charge is not greatly 
affected by natural convective fluid flow during 
typical vertical Bridgman growth of semiconductors 
(i.e,f materials with a small Prandtl number) [14rl51. 

(2) Transients. 

All models are quasi-steady; i.e.r the transient terms 
of the heat balance equations are assumed negligible. 

This assumption is normally valid due to the small 
lowering rates used in typical Bridgman growth. 
(Appendix E provides a brief discussion of teniperature 
transients.) 
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(3) Thermal coupling between the furnace and the 
charge. 

Thermal cabling between the charge and the furnace is 
described in the present work by a heat transfer 
coefficient, h. The heat transfer coefficient varies 
with position due, for example, to temperature 
dependent radiation heat transfer. The present thermal 
models assume that the heat transfer coefficient is 
constant within the hot furnace, cold furnace and 
gradient control region. Since the expected maximum' 
spatial variation of h within a furnace zone is 
approximately a factor of 2, it is considered that the 
variation of h is of secondary importance compared to 
the overall thermal coupling described by a suitable 
average value. Additionally, the hot and cold zone 
furnace temperatures are assumed to be uniform 
reflecting the heat pipe action. The lengths ot the 
furnace zones are assumed infinite. (Appendix D 
develops a criterion for the validity of this 
assunption) . 

(4) Constant thermal properties. 

Except for the change in thermal conductivity at the 
growth interface, all thermal properties are assumed 
independent of temferature. 
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2^ LITERATURE SURVEY 

t. number of thermal analyses of Bridgman-type crystal growth systems 
have been reported in the literature in recent years [9-11,17-28]. They 
can be classified according to the type of analysis which was performed, as 
shown belcw in Table 3.1: 


TYPE OF THERMAL ANALYSIS REFERE3X2S 


One-Dimensional (axial tem- 
perature variations only) 

9-11, 17-23 

Two-i)imensional (axial and 
radial temperature variations) 

17, 24-28 

Dimensional Variables 

10, 11, 18-20, 25, 27, 28 

Nondimensional Variables 

9, 17, 21-24, 26 

Analytical Results 

17, 22, 23 

Numerical Results (i.e., 
computer generated results) 

9-11, 18-a, 24-28 


Table 3,1; Current literature containing thermal analysis 
of Bridgman-type crystal growth systems. 


U. QNE-OIMENS T QNAL ANALYSES 

All of the one-dimensional m.odels describe the axial temperature 
distribution of a solidifying charge by equations similar to those used by 
Davis [18]. In these models, the charge is analogous to a fin moving 
through spatially varying boundary conditions. The equations for a moving 
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fin are found in many heat transfer texts [e.g.f 29,30]. The principal 
results obtained through a one*^imensional analysis are the axial 
temperature gradi«it at the growth interface and the axial position of the 
growth interface.. 

ii conjuction with his model, Davis [18] obtained experimoital results 
by solidifying load-tin alloys at various growth rates and hot and cold 
zone furnace temperatures. Good correlacion between experimental and 
calculated results demonstrated the usefulness of one-dimensional models 
for predicting the axial temperature distribution of the charge. Further 
experimental verificaticxi of one-dimensional modeling has been provided by 
Clyne [10,11] who solidified aluminum at various lowering rates. 

llie differences among the one-dim&isional analyses listed in Table 3.1 
concern several factors. The most apparent distinction is the type of 
system which is modeled. Some works deal with the furnace configurations 
of experimental growth systems employed for the growth of specific metals 
or metal alloys (Pb-Sn [18], A1 [10,19], Ag-Si [20]); others treat 
idealized furnace boundary conditions without reference to a specific 
charge material [9,17,21,22,231. The former analyses, being limited to 
specific charge materials, primarily analyze the effects of furnace and 
operation parameters such as the thermal coupling between the charge and 
furnace, position of the charge within the furnace, and lowering rate. The 
latter models are nondimensional and their results are more readily applied 
to an extended variety of growth systems. 

The one-dimensional analyses differ also with regard to the solution 
technique. Several employ finite difference or finite element methods 
[10,11,19-21]; others carry out analytical solutions which are sufficiently 
complicated to require computer evaluation of the resulting expressions 
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[9,18,22f23]. Only Chemg and Wilcox [17] and Sukanek [23] have been able 
to prcs«it results in simple analytical form. Generally^ the shift from 
numerical to analytical methods occurs as the modeling of the growth system 
becomes more idealized. 

PinallYf differences among the one dimensional analyses also pertain 
to the specific effects incorporated into the thermal models: inclusion of 
the crucible [10rll,18f20] r generation of latent heat at the growth 
interface [10,11,17,28,21-23] r unequal melt and crystal thermal 
conductivities [9-11,18-23], gradient control region [10,11,18,22,23,261, 
lowering rate of the charge [10,11,17,18,21-23], and temperature transients 
[ 10 , 11 , 21 ]. 

The majority of the one-dimensional models are concerned with the 
relationship between the position of the charge within the furnace and 
either the axial temperature gradient at the growth interface or the axial 
position of the growth interface. Chang and Wilcox [17] demonstrate that 
the interface position changes as the end of the charge approaches the 
region between the hot and cold zcne furnaces. Iheir results suggest that 
larger Biot numbers (see Section 5.2) and longer charge lengths tend to 
maintain a stationary interface location. Riquet and Durand [9] extend the 
results of Oiang ana Wilcox [17]. They show the manner in which the growth 
and lowering rates differ during a solidification experiment and indicate 
how the furnace temperatures should be varied with time in order to keep 
the interface location stationary. Sukanek [23] offers a simple criterion 
under which growth and lowering races will be nearly equal during 
solidification of a significant portion of the charge (see Section 5.5). 
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Two-dimensional thermal models are employed in order to obtain the 
interface shape. The difference between the one- and two-dimensional 
models is the incorporation of radial temperature variations. This 
addition^ however ^ makes the solution significantly more complex. Chang 
and Wilcox [17] used a highly idealized system (no crucibler gradient zone^ 
or geieration of latent heat> and equal thermal conductivites of crystal 
and melt) and obtained a solution in the form of two Fourier-Bessel 
infinite series. Evaluation of these series was carried out by computer. 
Sen and Wilcox [24], Fu and Wilcox [26], El-Mahallaway and Farag [25] and 
Domanski et al [27] used finite difference methods. Jones et al [28] used 
an electrical analogue model whose solution was also obtained by computer. 

A series of outstanding publications by Wilcox et al [17,24,26] 
demonstrate the usefulness of nondimensional variables in order to describe 
the thermal characteristics of a Bridgman growth system. Their results 
present two- dimensiaial isotherms in a solidifying charge for different 
values of several of these nondimensional variables. The effects of 
thermal coupling between the charge and the furnace (represented by the 
Biot number) and the lowering rate (Peclet number) was thus determined by 
Chang and Wilcox [17]. A significant conclusion is that the shape of the 
solidification isotherm depends strongly on its axial location within the 
furnace. Although not including a gradient control region in their 
thermal model, Chang and Wilcox [17] suggested the use of an insulating 
gradient control region in order to reduce isotherm curvature near the 
interface. Fu and Wilcox [26] later demonstrated the validity of this 
concept for a system which did not consider a crucible, latent heat, and 
conductivity change at the growth interface. Sen and Wilcox [24] showed 
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the effect of crucible conductivity on the isotherm shape for various Biot 
numbers; they also found that the thermal conductivity of the crucible 
should be close to that of the charge in order to maintain large axial 
temperature gradieits in the gradient zone. 

El-Mahallawy and Farag C25] calculated the two-dimensional temperature 
distribution in a solidifying Al-Al 3 Ni eutectic grown under various 
conditionsr including an abrupt change in the inner diameter of the 
crucible along the axial direction. Several comparisons of calculated 
results with experim^tal data showed good agreement 

lU EVALDATIQN OF 15E LITERATtlRE 

The cited literature provides a valuable background for further 
thermal analysis of Bridgman-type growth systems. In particular, the 
nondimensional formulation of Wilcox et al [17,24,26] identifies the 
significant combinations of thermal properties and system parameters which 
govern the thermal behavior of the system. Taken in part or as a whole, 
however, the previous work has not yet led to a systematic approach for 
Bridgman furnace design. The areas found deficient are; 

(1) Several important effects have not been studied, have been 
analyzed only in simple systems or have been addressed only 
for systems with specific thermal property values. For 
example, a general treatment of the behavior of radial 
temperature gradients near the growth interface is not 
available; only the case of a perfectly insulating gradient 
zone has received attention; generation of latent heat has 
been includ^jd only in one-dimensional models without 
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crucibles; the effects of unequal thermal conductivities of 
the crystal and the melt, and crucible thickness have been 
considered only for systems of specific property values . It 
is difficult to draw general conclusions regarding these 
effects. 

(2) Conclusions drawn from those analyses which did not use a 
nondimens ional formulatioi are difficult to exteid to systems 
of different parameter vadues. 

(3) As seen in Table 3.1/ all results, excepting parts of 
[17/22/231/ have been determined by computation. The 
presentation of these results is necessarily in graphical 
form and their interpretation quickly becomes awkward when 
several effects are acting simultaneously. In this respect, 
results presented in simple analytical form are more 
desirable. 
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lx THE Q^^E-DlMBSlQKftL lagESi&L MODEL 

The axial temperature distributicn of a solidifying charge determines 
two i jrtant characteristics of a particular growth experiment: 

Q) the axial temperature gradient in the melt at the growth 
interface. This must be greater than a critical value in 
order to prevent const’’ tutional supercooling (c.f., 
section 2.1.1) . 

(2) the axial location -'f the growth interface. In addition^to 
strongly influencing the interface shape, the temporal 
sequence of interface locations indicates the growth rate (as 
opposed to the lowering rate). 

Since neither of these factors pertains to radial temperature variations 
within the charge, the common approach has been to determine the axial 
temperature distribution through a one-dimensional thermal model, i.e., 
neglecting tne radial temperature variations. Such models neglect the 
radial thermal resistance within the charge compared to the thermal 
resistance between the charge and the furnace. The ratio of internal to 
external thermal resistance is expressed by the Biot number, Bi. A one- 
dimensional model is therefore valid for sufficiently small Bi (c.f., 
section4.2.3). 

The one-dimensional thermal model developed in this chapter posesses 
several unique features. The crucible and the radial thermal resisitance 
of the charge are accounted for in a way which readily demonstrates their 
effects on the axial temperature distribution without compromising the 
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simplicity afforded by the one-dimens icrial modeling. Additionally, small 
lowering rates common in semiconductor crystal growth is recognized to lead 
to simple analytical results from the one-dimensional model. Such results 
are used in Chapter 5 to demonstrate the effects of the various parameters 
of tile Bridgman growth system on the axial temperature distribution of the 
charge. 

AJ. THE MQVIMG FIN AgPBQACa 

A heat balance performed on a slice of the charge of infinitesmal 
length in the axial direction results in eq. [4.1] (see section 4.2.1) for 
the axial temperature distribution of the charge. The heat balance points 
out that the amount of heat conducted through the charge in the axial 
direction, and hence the axial temperature gradient in the charge, dep&ids 
directly on the amount of heat transferred to the charge in the hot furnace 
and from the charge in the cold furnace. Factors which affect the radial 
transfer of heat to or from the surface of the charge are therefore 
expected to affect the axial temperature distribution of the charge. 

The one-dimensional thermal model is initially developed in 
section 4.2.1 neglecting the crucible and the radial thermal resistance 
within the charge. The charge is, in this case, analogous to a fin moving 
through the spatially varying boundary conditions imposed by the furnace. 
The crucible, however, may provide a significant radial thermal resistance 
between the surface of the charge and the furnace boundary conditions. 
Further, a large Biot number indicates that the internal radial thermal 
resistance of the charge affects the radial heat exchange with the charge. 
These effects are approximated in sections 4.2.2 and 4.2.3 where it is 
shown that they can be included in the moving fin model as a modification 
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of the thermal coupling between the charge and the furnace (i.e., as a 
modification of the Biot number). 

THE DESCRIBING EQUATIONS FOR THE MOVING FIN MODEL 

The factors of concern for the development of a one-dimensional heat 
transfer model of the modified Bridgman growth system depicted in Fig. 2.2 
are shown in Fig. 4.1. The hot and cold heat pipes comprise the hot and 
cold zones; the region between them is called the gradient zone. The 
length of the charge is broken down into Lg, Lq, and L(; within the hot, 
gradient and cold zones respectively. The charge is lowered through the 
furnace with a velocity V, has crystal and melt portions with different 
thermal ccnductivities, and has a crystal-melt interface which generates 
lateit heat. A crucible provides containment for the charge. 

The thermal mode- makes the assumptions described in section 2.2.2. 
Additionally, the gradient zone is assumed adiabatic (i.e., h^ ® 0) in the 
one-dimensional model; this assumption is relaxed for the two-dimensional 
model of Chapter 6. 

Am2jl1 Chacge 

With the assumption of negligible radial thermal resistance within the 
charge, and neglecting the crucible, the charge is analogous to the moving 
thin rod treated by Carslaw and Jaeger [29] ; the equation describing the 
axial temperature distribution is; 

A 1^ = 0 14.1) 

where; A * cross sectional area of the charge 

P = surface area per unit length of the charge 
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where: Z > axial coordinate measured from the center of the 

gradient zone 

■ thermal conductivity 

j’ ■ mass density 

Cp ■ specific heat 

t ■ time 

Subscript "loc" serves to distinguish boundary 
conditions or charge properties of a particular 
furnace zone. 


T;n denotes the area**averaged charge temperature as a function of Z. The 
first term of equation [4.1] represents axial ccxiduciion of heat within the 
charge and the second term represents axial convection of heat due to 
motion of the charge at velocity V. The factor q" in the third term 
accounts for radial heat transfer to the surface of the charge^ per unit 
surface area of the charge. The fiml term represents heat storage due to 
temperature transients at fixed axial location^ Z. According to assumption 
3 in section 2.2.2, q" is given by: 

where: h » heat transfer coefficient between the furnace and the 

charge 

Tf ■ furnace temperature 


Substituting eq. [4.2] into eq. [4.1], neglecting the transient term 
according to assumption 2 in section 2.2.2, and nondimensionalizing the 
resulting expression, yields; 


j p d9iw 




[4.3a] 
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An equivalent nondimensional form of eq. [4.3a] is; 


j p ^ ^ Sf»\ 

d(Cf ■ ®i<.c iC.' 



[4.3bJ 


where! C'^CUi^) 

Equation [4.3b] has one fewer parameter than eq. [4.3aJ but the axial 
coordinater C / is less physical since the Biot number is included in its 
definition. 


4.2.2 The Crucible 

In Bridgman growth a charge is confined in a crucible which^ depending 
on conditions, varies in dimension and composition. Containment of the 
charge typically tends to decrease axial temperature gradients (Sen and 
Wilcox [24]). A crucible of low thermal conductivity lowers the axial 
gradient by adding thermal resistance between the charge and furnace, 
thereby decreasing the thermal coupling between the charge and the furnace; 
one of high thermal conductivity lowers these gradients by conducting heat 
transferred from the furnace axially within the crucible rather than within 
the charge. 

The domain of eq. [4.3] is considered to be the charge only and the 
crucible is not explicitly included. However, a model of the heat transfer 
within the crucible approximates its effect on q" of eq. [4.1] which is 
then expressed through modified Biot and Peclet numbers. 

Several previous one-dimensional models [18,22,23] have included the 
crucible only as a radial resistance between the furnace and the charge. 
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The present model for the crucible includes axial conduction as well as 
radial resistance by solving in an approximate manner the heat conduction 
equation for the temperature distribution of the crucible: 





0 


[4.4] 


where: f ■ radial coordinate 

^ ■ thermal diffusivity, Cp 

Subscript "cr" denotes the crucible. 

The first term of eq. [4.4] accounts for the radial thermal resistance of 
the crucible, while the second term accounts for the axial conduction of 
heat within the crucible. 

The principal assumptions involved in utilizing eq. [4.4] are that, at 
each axial location of the crucible, and 3']^^/Pf*are independent of 

the radial coordinate r and equal to d\,/62 and of the charge. 

(See section 4.4.1 for a discussion of these approximations.) The first 
term of eq. [4.4] can then be integrated to yield a radial distribution of 
Ter terms of the axial gradients in the charge. The boundary conditions 
at che crucible surfaces are: 




(at the outer [4.5a] 

crucible surface) 




(at the inner [4.5b] 

crucible surface) 


The heat transfer coefficient in eq. [4.5a] is between the furnace and the 
outer crucible surface. Equation [4.5b] assumes that the crucible and the 
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charge are in contact. The radial heat flow to the charge in eq. [4.1], 
q", is related to the radial distribution of by; 




3>Tc.r 


(at the inner [4.61 

crucible surface) 


The heat flow given by eq. (4.61 is substituted for q" in the third term of 
eq. [4.11. When nondimensiaialized, the resulting expression has the same 
form as eq. [4J] if Bi and Pe are replaced by their respective "effective" 
values, Bi* aix3 Pe*; 





\ 





' loc 


1 * 


[4.7] 


[4.8] 


where; % “ Dcr/D 

Kloc “ i^r/^^loc 

^^eC,l0C ® 

The relationship between the effective Biot number and the 
conductivity ratio K is shown in Fig. 4.2 for S * 1.25 and various Bi. 
It is seen that Bi* is significantly reduced by both low and high values of 
K, especially for larger Bi. In section 5.2 it is shown that increasing 
Bi* tends to increase the axial gradient at the growth interface. The 
conductivity ratio, K'*’, which maximizes Bi* for given S and Bi can be 
obtained from eq. [4.71: 
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♦ 



2 U "-0 




[4.9] 


For typical values of Bi and § , the conductivity ratio providing the 
maximum effective Biot number is found to lie between 0.1 and 1.0. 

It is sometimes possible to increase Bi* by increasing the thickness 
of the crucible (i.e., by increasing S ). This effect is similar to the 
"critical radius of insulation" effect whereby a wire which is conducting 
electricity may be able to dissipate more heat to its surroundings with a 

i 

layer of insulation on its surface [301. For the charge/crucible system, 
an increase in S increases the charge surface area expose -.o cat 
transfer with the furnace. If this effect more than offsets the o^si.ig 
effects of increased radial resistance and axial conductance of the 
crucible, then Bi* increases with S . 

Consider that Bi* in eq. [4,7] is a function of S with constant K 
and Bi. The value of S which maximizes Bi*, S'*" , is given by tho roots of 
the following equation: 


Equation [4.10] has positive rotts greater than unity ( S > 1 when a 
crucible is present) only when K < 1/2 and Bi < 1/4. Figure 4.3 plots the 
values of h* for various K and Bi. Tt is seen that it is advantageous to 
have a thicker crucible when Bi is small, but only if K is also small. 
Figure 4.4 shows the increase in Bi* when S “ S'* compared to a bare 
charge (i.e., compared to S ■ 1). (It is found that this ratio is not a 
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Strong function of K but correlates well with S"^.) Accordingly, for those 
cases where an increase in S results in an increase in Bi*, the maxiirum 
effect is not large. Therefore, increasing the crucible thickness should 
not be considered as a viable method for increasing the axial temperature 
gradients near the growth interface. 

i.2^ Radial Temperature Gradients Wit hin the Charge 

Equation [4.1] considers that the temperature within the charge at 
each axial location is uniform. Presently, radial gradiaits in the charge 
are considered insofar as they affect the thermal coupling between the 
furnace and the charge and thereby the axial temperature distribution. 

When radial temperature variations in the charge are taken into 
consideration, eq. [4.3] assumes the form (see Appo:idix A): 

IZ? ' ^ 9") - ( 9^ - = 0 I4.ni 


where: » surface tenperature of the charge, i.e., at ^ * 1/2 

The term 93 - 9^^ in eq. [4.11] accounts for the effect of radial 
temperature variations within the charge on the q" term of eq. [4.1], An 
approximation for 9s - 9m is obtained by solving in an approximate manner 
the two-dimensional heat conduction equation in the charge: 


11 - 

f 






= 0 


[4.12] 


Assuming that the Peclet term is negligible (see section 5.1) , and, as in 
the crucible effect, that is only a weak function of the radial 
coordinate, ^ , eq. [4.12] can be integrated in the radial direction with 
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i*9w/(jC‘as a variable. The integration indicates a parabolic variation for 
the charge temperature at a particular axial location; 


- - 






and: 


9-0 = - — 


[4.13a] 


[4.13b] 


Combining eqs. [4.11] and [4.13b] leads to a new effective Biot number, 
Bi**, replacing Bi in the third term of eq. [4.3a], which accounts for 
radial temperature gradients within the charge as well as within the 
crucible; 



I 



[4.14] 


According to eq. [4.14], Bi** = Bi* for Bi* « 8; i.e., the radial 
thermal resistance within the charge does not affect the axial temperature 
distribution. Equation [4.14] indicates also that the maximum value that 
Bi** can attain is 8 when Bi* is very large. The internal radial 
resistance of the charge thus limits the degree to which axial temperature 
gradients can be increased by an augmentation of the thermal coupling 
between the charge and the furnace. 

The validity of eq. [4.14] may be tested by comparing the axial 
temperature distributions for a system whose internal thermal resistance in 
the radial direction is considered (1) exactly through the use of a two- 
dimensional model and (2) through the use of eq. [4-14]. A stationary fin 
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(without crucible) of infinite length provides a simple system for which 
analytical results for both models are readily available (see Fig. 4.5). 
In order to consider the most extreme case, the heat transfer coefficient 
at the surface of the fin is taken to be infinite; for these conditions, 
Bi is infinite and Bi** * 8. in this case, the expression for the cross 
sectional area-averaged temperature for the two-dimensional fin is (31): 



n»i 


where: 

= 0 

J*q » Bessel function of order zero 

and the axial temperature distribution for the one-dimensional fin is [31] : 

[4.16] 

Equation [4.16] with Bi** » 8 is compared to eq. [4.15] in Fig. 4,6 with 
good results. Equation [4.14] is therefore taken to be a reasonable 
anaytical approximation for the effect of radial gradients in the charge on 
the overall axial temperature distribution. 

Abrupt changes in surface boundary conditions affect the temperature 
of the charge more at its surface than at its center. The assumption of 
constant 3^9/SC^ over the cross section of the charge and, as a 
consequence, eq. [4,14], are expected to be less valid near such changes in 
boundary conditions. The moving fin model considers an abrupt change in 
the heat transfer coefficient at the ends of the gradient zone. Equation 
[4.14] is least accurate near these locations but becomes increasingly 
valid with distance into the furnace zones. A further verification of 
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9t = 0, Bi**oo 




(b) One-dimensioncil model used with eq.[4.14] 




c 


Fia. 4.5: >iodels of a simple stationcur..' fin 'jsed to check the validity 

of eq. [4.14] . 
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eq. [4.14] obtained for the conditions of an adiabatic gradient zone is 
provided in section 4.6. 

somnoN OF moving fin squatiqns 

Using the results of sections 4.2.2 amd 4.2.3 in order to include the 
effects of the crucible and radial thermal resistance within the charger 
eq. [4Ja] becomes: 

The parameters Pe*r Bi**, and ©f are not constant throughout the furnace. 
According to assumptions 3 and 4 in section 2.2.2, however, they are 
constant within each of four separate furnace regions: the hot and cold 
zones and the liquid and solid parts of the gradient zone. (It is assumed 
that the growth interface is located within the gradient zone.) The 
solution of eq. [4.17] within each region yields two exponential terms for 
the homogeneous solution and a particular solution that depends on ©f. The 
assumption of uniform hot and cold zone furnace temperatures yields a 
simple constant for the particular solution in tiiese zones. A particular 
solution is not needed in the gradient zone since hQ = 0. 

Hie constants of integration for the homogeneous solution of eq. [4.17] 
are found using the boundary conditions of equality of temperature and 
continuity of heat flux between adjacent regions. The temperature boundary 
condition is: 

where the subscripts "j" and "k" refer to opposite sides of a boundary 
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between adjacent regions. The flux boundary conditions Uoed in the moving 
fin model are: 


<j9»n,W 

AC, ^ AC, 


between the hot and 
gradient zones 

— or — [4.19] 

between the gradient 
and cold zones 


[Rk*K.U--0](^\ = 

Section 4.4J explains that the axial gradients in the gradient zone which 
result from the use of eg. [4.201 should be interpreted as average 
charge/crucible axial gradients; these average gradients will be denoted by 
Gl and Gg for the liquid and solid sides of the interface respectively. 
Unavoidable radial heat transfer between the charge and the crucible near 
the growth interface causes Gl and Gg to significantly differ from their 
respective counterparts in the charge, and . The 
relevant axial gradient for the supercooling requirement is, however, that 
in the charge. Evaluation of the difference between G and at the 
growth interface requires two-dimensional modeling and is therefore 
considered in Chapter 7. 

Boundary conditions at the end of the charge must be applied when the 
charge length can not be considered infinite (c.f., section 5.5). 
Appendix D discusses the calculation of approximate Biol '•bers for the 
nondimensional end boundary conditions. 

The constants of integration have been determined analytically for 
systems of infinite charge length; expressic's for the axial gradient in 





at the growth 
interface 


[4.20] 


- 63 - 


ORlGlNAt PACE B 
OF POOR QUAUTY 


the melt at the interface and for the axial position of the interface of 
such systems ere given in Appendix B. Of special interest are the 
simplified results obtained for the small lowering rates common to 
semiconductor crystal growth, eqs. [B.9] and [B.IO]. If the complete axial 
temperature distribution is desired, or if the charge length is not 
infinite, it is more caivenient to determine the constants of integration 
by computer. 


AA ASSUMPTIONS IN THE 




3 fIM MODEL 


The validity of the primary assumptions employ^ in the development of 
the moving fin thermal model of section 4.2 is discussed in this section. 


4»4.I Equal Axial Gradients of the Charge and Crucible 

In order to obtain a simple analytical expression for the effect of 
the crucible thickness and conductivity on the axial temperature 
distribution of the charge, the crucible model (section 4.2.2) assumes that 
the first and second axial temperature derivatives in the crucible are 
independent of the radial coordinate, , and equal to their re>ipective 
values within the charge. These approximations are motivated by the two- 
dimensional temperature distributions presented by Wilcox et.al. [24,26]. 
Sen and Wilcox [241 determined isotherm shapes in a charge/crucible system 
for various K and Bi by a two-dimensional finite difference computer model 
under the following simplified conditions; ^ * 1*25, 

.75, Rj( * 1.0, Rh = 0.0. The shapes of their numerically 
computed isotherms indicate that the assumptions are reasonable for values 
of K greater than about 1/2. As K decreases below 1/2 (e.g., as k^r 
decreases), increasing radial resistance of the crucible produces a 
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pronounced radial temperature variation in the crucible; the axial 
gradients in the crucible are larger than those in the charge near the ends 
of the hot and cold zone. It is therefore expected that the crucible 
effect, as described by eqs. [4.7] and [4.8] will be most in error for 
small values of R. 

Ihe temperature distributions of Fu and Wilcox [26] indicate that the 
radial variation of the axial gradients are most pronounced at the ends of 
the hot and cold zones, especially for systems without a gradient zone. 
The absence of a gradient zone in the model of Sen and Wilcox [24] 
therefore accentuates the error in the approximation for systems with a 
gradient zone. 

The two-dimensional thermal model developed in Chapter 6 does not 
assume that axial temperature gradients in the charge and crucible are 
equal. The effect of the crucible on the axial temperature gradients of 
the ch£urge determined by the two-dimensional model are compared to results 
given by the one-dimensional model using Bi**. Agreement is good even for 
values of K as small as 0.1. It is therefore concluded that the effective 
Biot number provides a reasonable estimate of the crucible effect on the 
axial temperature distribution of the charge. 

AAm2 Equal Growth and Lowering Rates 

The term PegRH ^ [4.201 represents the generation of latent heat 
at the growth interface; the velocity used in Pes should therefore be the 
actual growth rate, R. When the charge length is not infinite, however, the 
growth rate exceeds the lowering rate, V, [9,11,221. The present model 
neglects the difference between V and R under such conditions and uses Pes 
based on the lowering rate V. in this way, the mathematical solution is 
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simplified since the rate of genera tic« of latent heat then has a known and 
constant value. On the other hand, if the growth ratCf Rr were used, the 
rate of generation of latent heat would not be known a priori or constant. 

The assumption that R«V in not exessively restrictive. Riquet and 
Durand [91 and Sukanek 122] show that the growth rate is most different 
from the lowering rate when the ends of the charge are passing zone 
boundaries. When the length of the charge in both the hot and cold zones 
can be considered infinite (see section 5.5), the growth rate and lowering 
rate are equal. Further, growth conditions may be such so that the effect 
of latent heat on the axial temperature profile is negligible (see 
section 5.4); the error created by the use of V rather than R in Pes would 
then be likewise small. Lastly, if desired, the solution of eq. [4.171 can 
be modified so that Pes in eq. [4.201 is based on R as follows: The axial 
position of the growth interface for two closely spaced locations of the 
charge in the furnace is calculated with an initial assumed value of R. 
The difference in position of the interface for these two locations 
indicates a growth rate which need not be equal to the original assumed 
value. The growth rate is accordingly modified and the procedure iterated 
until convergence is achieved. 

4.4.3 Flux Boundary Condition at the Growth Interface 

Consider that the charge and crucible temperatures, and 
respectively, are each represented a single value at each axial 
location. (The variable 4> is used in this work to denote non-dimensional 
temperatures employed in two-dimensional thermal models.) When PesRn ^ 0 
and/or Rr C, the axial gradient in the charge, r must undergo an 

abrupt change at the growth interface. The axial gradient in the crucible. 
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on the other hand, does not change at the growth interface. Figure 4,7 
shows that such conditions lead to a difference between the charge and 
crucible temperatures as well as their axial temperature gradients. The 
temperature difference leads to an exchange of heat between the diarge and 
crucible near the growth interface. The relationship between this radial 
heat transfer and the shape of isotherms is discussed in section 7.3; the 
pres&it section addresses its effect on the flux boundary condition at the 
growth interface for the moving fin model, eq. [4.20]. 

Define an average charge/crucible temperature weighted by their 
respective axial thermal conductances: 


^ [4.21] 

kioc is respectively and ks cor the liquid and solid portions of the 
charge. Nondimensionalizing eq. [4.21J yields: 




[4.221 


Differentiating eq. [4.22] in the liquid and solid yields: 



[4.23a] 




[4.23b] 


• Ml AlKl ftt _ page 



Fig. 4.7 ; Crucible and mean charge terperature distributions 
near the interface vvhen Pes^H ^ and/or Rj^ # 1. 
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The moving fin model assumes that E q " 0. In th..s casSf theie is no heat 
transfer to or from the surface uf the crucible in the gradient zone and G 
must be invariant in the gradient zone except for a step change at the 
growth interface. The relationship between Gl and Gs is found by 
considering the heat flux boundary condition at the growth interface. In 
nondimensional form: 


Rk 


c)<^ \ _ 

li;i “Ua 




[4.24a] 







, [4.24b] 


where subscript "i" denotes the interface 


By multiplying eg. [4.24b] by Ks(S*-i) ^ adding to eg. [4.24a] and using 
egs. [4.23] r the following relation is obtained: 

[Rk* K5(S^-0]Gi=[|-*KsU'-''>]G, * 14.251 


The one^imensional moving fin model neglects the difference between 
and / i.e.f this caser egs. [4.23] show that: 



[4.26a] 


G 


3 

s 



[4.26b] 


Substituting egs. [4.26] into eg. [4.25] results in the interface boundary 
condition used by the moving fin model, eg. [4.20]. These results suggest 
that the axial temperature gradient in the gradient zone predicted by the 
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moving fin models t should be interpreted as the average 
charge/crucible axial gradient, G. The difference between G and is 
discussed in section 7.3. 

4.4.4 Diabatic Gra dient Control Region 

The moving fin model assumes that there is no heat exchange betw‘ en 
the charge and the furnace within the gradient control region; i.e., 
h (3 ■ 0. Chang and Wilcox [17] suggested the use of an adiabatic gradient 
control region as a means of reducing isotherm curvature near the growth 
interface. Further motivation for this assumption is that it sijnplifies 
the analytical expression of the temperature distribution in the 
solidifying charge. 

A perfectly insulating gradient control region is a conceptual limit 
and is unobtainable in practice. Radial heat transfer to or from the 
charge in the gradient :-.one affects isotherm curvature; the assumption of 
an adiabatic gradient control region is therefore eliminated in the two- 
dimensional model described in Chapter 6. The present section formulates a 
simple model of heat exchange within the gradient control region in order 
tc determine when this heat exchange is negligible with regard to the axial 

temperature distribution of the charge. 

Ccxisider the portion of the fumace in the gradient control region to 
be an annulus located between and in contact with the hot and cold heat 
pipes. The gradient zone annulus has conductivity cross sectional area 
Aq, and has a temperature, TqCZ) which depends only on the axial 
coordinate. Heat exchange between the gradient zone annulus is described 
by the gradient zone heat transfer coefficient h^. No heat loss is assumed 
at the outer diameter of the gradient zone annulus. (This model of the 
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gradient control region is identical to that used in Chapter 6; c.f.# 
section 6.2 and Figs. 6.1 and 6.2.) 

In order to develop a relatively simple expression for the temperature 
distributions of the charge and gradient zone annulus, the present system 
is specified by the following parameter values: 

1. Equal heat transfer coefficients in the hot and 
cold zones. 

2. Axial convection of heat due to the charge motion 
is negligible; i.e., Pe is small. (See 
section 5.1.) 

3. The effect of the generation of latent heat is 

negligible? (Pes% “ section 5.4.) 

4. Equal thermal conductivities of the melt and 
crystal ■ 1). With condition 1 above, 

_ . ** _ . h* 

BiH - Blc . 

5. Equal charge lengths in the hot and cold zones 

) • 

A system which has the above parameter values is referred to in the present 
work as a "symmetric" system. The axial temperature distribution of the 
charge for a synmetric s^'stem is: 

It is only necessary to consider one-half of a symmetric system, 
either > 0 or C, < 0. An additional assumption used in this section 
is that the charge length is infinite; in this way, the absence of the end 
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boundary coiditions and charge length as thermal parameters simplifies the 
expressions for the temperature distributions of the charge and gradient 
zone annulus. 

A schematic of the temperature distributions in the charge and 
gradient zone annulus is shown in Fig. 4.8. Only the cold portion of the 
symmetric system is shown. At the center of the gradient zone (i.e.r at 
« 0) , the assumpticxi of a symmetric system dictates that ©m * 9 g “ 

The cold end of the gradient zone annulus is at the cold zone heat pipe 
tai^rature. 

The equation governing the mean charge temperature in the cold zcxie is 
adapted from eq. [4.17] assuming that Pe = 0.0: 

^ ^ [4.28] 

where: ^ ® 

The charge temperature in the gradient zone is described by eq. [4.28] with 
Bic* replaced by BIq* and 0f^c replaced by 

i.” (.64 - > 0 14.291 

where: 0 ^ ^ ^ ^ 

The fin equation for the gradient zone annulus can be derived by noting 
that the transfer of heat between it and the charge/crucible requires the 
following relation between the second axial derivatives: 
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Fia 4.8; Scheitiatic of the temperature distribution in ^e charge and 
^adiSA zone annulus for a syitmetric system '.dien Big 0. 
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(Equation [4.301 assumes that the axial conduction of heat in the crucible 
is described by the axial derivatives of 9^ .) Substitucinq eq. [4.301 
into eq. [4.291 : 


AC.’- 







where; 


KslS'-'') = 



[4.311 


A 

Equations [4.281, [4.291 and [4.311 are simultaneously solved 
utilizing th<s i.wllowing boundary conditions: 


e^io) -o.s 

0 


= 0.5 

0W, Cop') - 0 


[4.321 




L 




• <h 

L ac 1 


The resulting distribution for the temperature difference between the 
charge and gradient zone annulus is; 
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K 


R|_ 


U.33J 






■tank (^Wj) 




where: = 2 [BI** (\ + 

V. = 2 (Blf 

^4 = R.i/D-'KrU'-n] 

k 

Equation [4e33] demonstrates that the difference in temperature between the 
gradient zone annulus and the charge varies exponentionally and decreases 
more rapidly for larger BIq* and smaller Rq# A large value of Big* is 
indicative of effective thermal coupling between the charge and gradient 
zone annulus, A small value of Rq indicates that the gradient zone annulus 
has a relatively small capacity to conduct heat in the axial direction 
compared to that of the charge/crucible combination; consequently^ any heat 
transfer to or from the gradient zone annulus serves to rapidly equilibrate 
®<i, to . 

Uie importance of heat transfer in the gradient zone compared to that 
orcuring in the hot and cold zones on the axial temperature distribution of 
the charge is revealed by the expression tor the axial temperature gradient 
in the charge at the center of the gradient zone; 




^ 1 

i 


1 1 "1 

\ w 

RM 1 

^<9 


1 1 


* ^ M 

'*■ 1 tink ^ 


[4.34] 
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Figure 4.9 shows the variation of G( C * 0) with systems with 
Bi^* * Bic*. When Rq = C, the gradient zone does not conduct heat from the 
hot or cold zones for subsequent excnange to or from the charge; in this 
case, the gradieit zone may be called adiabatic even though Bi <3 0. As Rq 
increases, the temperature distribution of the gradient zone annulus 
becomes increasingly linear providing a temperature difference between it 
and the charge; the resulting heat exchange increases the axial temperature 
gradient of the charge. As Rq approaches infinity, the heat exchange 
between the gradient zone annulus and charge is little affected by 
changes in Rq since 6 g( C ) approaches its limiting linear variation. 

Figure 4.10 demonstrates the effect of BIq* for systems with 
Rq = 1.0. The adiabatic gradient zone case is indicated by the curve 
labeled BiQ*/Bic* » 0.0, Heat exchange between the charge and gradient 
zone annulus becomes increasingly important as Big* increases relative to 
Bic*f resulting in a significant increase in the axial gradient of the 
charge. In the limit as Bic* approaches infinity, the temperature 
distribution of the charge in the gradient zone varies linearly between the 
boundary conditions imposed by the hot and cold zone furnace. 

The results presented in this section were based on a model which 
assumed infinite charge length. Further, Figs, 4.9 and 4.10 examined the 
special case of = 1.0. It is expected that the effect of a diabatic 
gradient zone increases as the relative length of the gradient zone 
increases (e.g., finite charge length and increasing Nevertheless, 
the results of this section suggest that the gradient zone can be made 
effectively adiabatic by two design options: 


(1) Small Rq. a material of poor thermal conductivity placed 
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between the hot and cold zone furnaces provides for a small 
value of R<;. 

(2) Small BIq relative to the Biot numbers in the hot and cold 
zones. This may be accomplished in systems where the 
conductive component to the gap heat transfer coefficient is 
small (see Appeidix Q by placing highly reflecting radiation 
shields between the gradient zone annulus and the charge. 

THERMAL PARAMETERS QF TBE MOVING FIN VCLZL 

The thermal parameters which appear in the moving fin model are 
sumnarized in Table 4.1 

AjA (3>gARI9QN WI TH THE TWP- DIMENSIQNAL RESULTS QF FU AND WILQQX [26] 

Results derived from the preset one-dimensional model can be compared 
to similar results of Fu and Wilcox (26] who included radial temperature 
gradients within the charge in a two-dimensional finite difference model. 
Figure 3 of their paper presents the radial variation of the axial gradient 
at the growth interface for systems of infinite charge length in the 
absence of a crucible. Point values, estimated from the curves of their 
Fig. 3, were numerically integrated to obtain the average axial gradient 
over the cross section at the growth interface. These results are compared 
to the values obtained using eqs. [4.14] and [5,8] for several values of 
Bi and in Table 4.2. 

Since the model of Fu and Wilcox [26] does not include a crucible, the 
comparison shown in Table 4.2 tests the validity of eq. [4.14] in 
accounting for the effect of radial thermal resistance of the charge on the 
effective Biot number. Considering the approximations required in the 
derivation of eq. [4.14], the agreement of the data must be taken as 
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excellent. It is considered, therefore, that eq. [4.14] provides a 
reasonable correction for the effect of radial temperature variations 
within the charge on the axial temperature distribution of the charge. 
(The validity of eq. [4.7] in accounting for the effect of the crucible cn 
the effective Biot number is provided in section 6.3.; 


parameter 

definition 

comnents 

Bifl 

hflO/kL 

" , 

Bic 

h<p/ks 


Pes 

VDkg/ f Cp 

PeL * Pes^K 

s 

Dcr/0 


KS 

ICcr/*^S 

Kl * Ks/% 

Ro<,s 



Roc,l = Ro.,s/% 

% 

ki/ks 1 

% 

AHsi/Cp(Tf,H “ Tf,c) 



I<3/D 

1 

j 

x» 

Ifl/D 

required only if the 
charge can not be 
considered infinitely 
long (c.f., section 4.2.5) 

Xc 

Lc/D 

^^H,end 

^H,en(p/^L 

®^C,end 

hc,en<p/ks 

1 

Cl 

Zi/D 


©i 

(Ti - Tf,c) 

only one is independent 


Table 4.1: Thermal parameters in the moving fin model. 
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CSO/SC'Jrtv 

G 


Fu and Wilcox [26] 

Eqs. [4.14] and [5.8] 


0.0 

1.310 

1.261 

Bi ■ 2 

0.25 

0.975 

0.961 

(Bi** - 1.6) 

0.50 

0.778 

0.775 

Bi - 0.4 

0.0 

0.616 

0.617 

(Bi** - 0.38) 

0.5 

0.464 

0.472 


Table 4.2: Axial temperature gradients at the growth interface, 
Q , as obtained through eqs. [4.14] and [5.8], and 
average axial temperature gradients, QS/5 C;')av , as 
determined from the two-dimensional model of Fu and 
Wilcox [26]. (Symmetric system and no crucible.) 
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1m, TSL axial temperature DTSTRIBtTTTQN O F THE CHARGE 


This chapter describes the effects of the thermal parameters of the 
vertical Bridgman growth system on the one-dimensional thermal behavior of 
the charge. The presentation of the results emphasizes the thermal 
reguirements for satisfactory single crystal growth developed in 
section 2.1 which relate to the axial temperature distribution in the 
charge: (1) the axial gradient in the liguid at the growth interface 
(constitutional supercooling) and (2) the axial position of thv* growth 
interface (a factor strongly affecting the shape of the interface). 

Axial temperature distributions in the charge are calculated bv 
computer using the moving fin model described in (2iapter 4. Such results 
indicate how the interface position (or the nondimensional melt 
temperature) is affected by t.he thermal parameters. The axial temperature 
gradient within the melt in the gradient zone, is studied by utilizing the 
analytical results presented in Appendix B. 

Section 4.4.3 explains that the axial temperature gradient at the 
interface calculated from the moving fin model must be interpreted as the 
average of the charge and crucible axial gradients. The variable G is used 
to denote this average value. The difference between G and the axial 
temperature gradient in the charge is discussed in section 7.3. 

For a particular growth experiment having constant process pa rameters 
(e.g., furnr.ce temperature, lowering rate) the temperature of the interface 
and its location within the furnace are functionally related. Since the 
shape of the interface strongly depends on its axial location 117,26 and 
Chapter 7], it is expedient to consider the interface location as the 
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independent variable. Given the interface location, the thermal model 
predicts a nondimensional interface temperature, ei; the dimensional 
interface temperacure, Tj^, is related to the hot and cold zone furnace 
temperatures and ©i as follows: 

The desired interface location is achieved by adjusting the hot and/or cold 
zone furnace temperatures so that eq. 15.1] is satisfied. 

In consideration of the many variables in the moving fin model, the 
synmetric system defined in section 4.4.4 is used as a reference against 
which the effects of individual parameterr will be compared. 




The Peclet number expresses the ratio of axial heat transfer in the 
charge due to the lowering velocity, V, and due to ronduction. The 
effective Peclet number, Pe*, incorporates the effect of the crucible on 
the convection term of eg. [4.1al. 'Hie unmodified Peclet number is used to 
describe the generation of latent heat (eq. [4.20]). 

The effect of Pe* on the axial temperature distribution in an 
otherwise symmetric system obtained through eq. [4.3,7] is shown in 
Fig. 5.1. It is seen that Pe* tends to increase the temperature of the 
charge at all locations; the effect is more pronounced as Pe* increases 
(e.g., high lowering rate or low tl,-. -nal conductivity of the charge). As a 
consequence, in systems with fixed hot and cold zone furnace temperatures 
(i.e., the nondimensional melting point temperature is constant), the 
growth interface moves toward the cold zone as Pe* increases; 
alternatively, the nondimensional melting point temperature must increase 
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(i.e., one or both of the furnace temperatures must be lowered) if the 
interface is to remain at the same axial location as Pe* increases. 

Table 5.1 gives expressions for the Peclet number based on solid 
properties for Ge and CdTe (materials with high and low thermal 
conductivity, respectively) as a function of lowering rate V and charge 
diameter D. Using lowering rates commcxi to Bridgman growth experiments 
(e.g., .1-10 ^ m/sec), it is seen that the Peclet number is small compared 
to unity. A test criterion for its relative magnitude is provided through 
the characteristic roots of eq. [4J.71: 


W. 






•+ 1 (f? 


Bi" ■ 


'li 


I 


[5.2] 


Expanding eq. [5.2]: 





1 1 

I 

2 

Z igBl** 8 


+ • • • 7 


[5.3] 




Peg 

Ge 

.097 cmVsec 

VD/967 

CdTe 

.011 cm2/sec 

VD/111 


(X = thermal diffusivity 
V measured in ^m/sec 
D measured in cm 

Table 5.1: Magnitude of typical Peclet numbers, 
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If the following inequality holds, the characteristic values, are of 
opposite sign and their sum is equal to Pe*: 


if: 


1 (,Pe*f 

2 


1 


then; 



2 


[5.4] 


[5.5] 


Further, Pe* disappears from the characteristic roots under the following 
conditions: 


if: 


then: 


LS « ] 

±z 


[5.6] 

[5.7] 


Equation [5.6] agrees with Chang and Wilcox [17] who reported that the 
Peclet effect is stronger for smaller Biot numbers. 

If the inequality in eq. [5.6] holds, the Peclet effect is small enough 
to satisfy the Peclet number criterion for symmetric systems (i.e., Pe=C); 
in addition, simplified analytical expressions are obtained for the axial 
gradient within the melt in the gradient zone (Appendix R), Since 
eq. [5.6] does hold for typical Bridgman growth, thf \mplified 
expressions are used in the remainder of this chaiJter. 

Equation [4.20] shows tlriat the gaieratioi of latent heat at che growth 
interface is also dependent on Pe. The charge lowering rate, therefore, 
has a small effect on the axial temperature distribution if both eq. [5.6] 
holds and the latent heat effect is small. (Sections 5.4 and 7.3 address 
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the effect of latent heat on the axial temperature distribution in the 
charge and on the radial temperature variations near the interface, 
respectively. ) 

BIOT NUMBER EFFECT 

l^e Biot number, through the heat transfer coefficient, h, is a direct 
measure of the thermal coupling between the charge and the furnace. The 
effective Biot number, Bi**, incorporates the effects of the crucible and 
radial temperature gradients within the charge on this thermal coupling. 
Typical values for the effective Biot number vary from 0.05 for high 
conductivity materials such as Ge to 5.0 for low conductivity materials 
such as CdTe. (Appendix C gives approximate values for typical heat 
transfer coefficients.) Axial temperature profiles for several Bi** are 
shown in Fig, 5.2, In agreemeit with the results of Chang and Wilcox [17], 
it is found that the charge temperature follows more closely the furnace 
temperature and, as a result, the axial temperature gradient in the 
gradient zone increases, as Bi'* increases. 

The expression for thfe axial gradient near the growth interface 
presented in Appendix B for infinite charge length can be simplified for 
symmetric systems; 

The dependence of the axial gradient on Bi**, according to eq. [5.8] , is 
plotted in Fig. 5.3 for various gradient zone lengths, . The curves 
show that the dependence of the axial gradient on becomes stronger with 
increasing Bi**. It can also be seen (curve = 0) that there exists a 
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mAlu 


minimum Bi** for any desired nondimensionai axial gradient. 

The Biot numbers of the hot and cold zones are generally not equal 
due, for example, to the temperature dependence of radiative heat transfer. 
For such conditions, the zone with the larger effective Biot number more 
strongly influences the overall temperature level (Fig. 5.4). Ctompared to 
the symmetric case where Big** » Bi^**, the charge temperature for 
asymmetric systems increases when Big** > Bi^** and decreases when 
Bic** > Big**. If the location of the crystal-melt interface for each 
curve of Fig. 5.4 is considered to be the center of the gradient zone, the 
effect of unequal Biot numbers appears as a change in nondimensionai 
interface temperature. To retain the interface in a given location 
requires lowering the hot and/or cold zone furnace temperatures as Big** 
increases relative to Bi<;**. 

The expression for the axial temperature gradient in the gradient 
zone for a system with infinite l^gth which is symmetric except hg ^ 
is; 



[5.9] 


Equation [5.9] indicates that the smaller of the effective Biot numbers has 
a greater effect on decreasing the axial gradient. Efforts to increase 
axial gradients by adjusting the heat transfer coefficients should 
therefore first be directed at the zone with the smaller Bi**. 
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^ EFFECT QF UN EQUAL THERMAL CONDUCTIVITY OF THE CRYSTAL AND TBE MELT 

Rk expresses the ratio of melt to crystal thermal conductivity. 
Semiconductors have Rr values greater than unity whereas, for metals, the 
value of Rr is less than unity. The effect of Rr on the axial temperature 
distribution in the charge obtained through eq. [4.17] is demonstrated in 
Fig. 5.5 for systems with S ■ 1 and which are otherwise symmetric (i.e.. 
Big ■ Bi<;/RR). Hie charge phase with the higher thermal conductivity tends 
in all instances to have a lower axial gradient because of lower thermal 
resistance to heat transfer in the axial direction.- Hie axiad gradient in 
the melt near the interface is therefore less than in the crystal for 
charges v- ^th Rr > 1. 

The expression for the axial gradient in the melt near the growth 
interface for a system of infinite length which is symmetric except Rr ^ 1 
is (^pieidix B): 





C5.10] 


A comparison of eqs. [5.4] and [5.10] shows that a value of % greater than 
unity is detrimental to the establishment of large axial gradients in the 
melt. The axial gradient as a function of Rr according to eq. [5.10] is 
plotted in Pig. 5.6 for several values of Bic. It is found, for example, 
that this conductivity effect may reduce the gradient in a germanium melt 
(Rr = 2.5) by about 50%. 
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The quantity of latent heat of solidification liberated at the 
crystal-melt interface is given in nondimensional form by the product Pes% 
(eq. [4.20]). The expression for the axial gradient in the melt near the 
growth interface in a system of infinite length which is symmetric except 
PcsRh 0 is: 


r _ O-Bi.**) Xa ~2Ci1 

It can be seen that the effect of latent heat on the axial gradient (and 
also on the axial temperature profile) is small if: 


Equation [5.12] demonstrates that the effect of latent heat on the axial 
temperature behavior is larger for smaller Bi** and for larger lowering 
rates of the charge (Peg is proportional to V). 

Axial temperature profiles for various values of ^ are given in 
Fig. 5.7. It is seen that the generation of latent heat both increases the 
charge temperature and decreases the axial gradient in the melt; the latent 
heat effect disappears for small values of . 

(The effect of latent heat and the Peclet number effect described in 
section 5.1 are coupled through their mutual dependence on the lowering 
rate, V. In order to isolate the latent heat effect in Fig. 5.7 , Pe is 
chosen sufficiently small so that eq. [5.6] is satisfied. Such a small- 
value for Pe would normally also eliminate the latent heat effect; i.e., 
eq. [5.12] would also be satisfied. The values of in Fig. 5.7 were 
obtained oy choosing appropriately large values of %.) 
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^ "INFINITE* CHARGE LENSm 

When the charge is infinitely long/ the temperature field of the 
charge does not channe during growth and, hencer the interface position 
remains fixed. Figure 5.8 shows the progression of axial temperature 
profiles as the charge is lowered from the hot zone to the cold zone for 
charge loigths less then infinite. It can be seen that charge temperatures 
are displaced toward the hot zone furnace temperature when most of the 
charge is in the hot zone and vice-versa. Accordingly, to achieve ccnstant 
interface position for non-infinite charge lengths, the nondimensional 
solidification temperature, must be reduced as the experiment proceeds. 
At ccnstant 9^, the growth rate is greatar than the lowering rate while the 
interface moves from the cold to the hot zone. 

■nie contribution from the positive characteristic root, eq. [5.2], to 
the solution of eq. [4.17] is normally small and is zero for an infinite 
charge length. The charge thus appears infinite in length when the 
contribution from the negative root also becomes small, i.e., for large 
Ct . The temperature change within the hot or cold zone reaches 
approximately 99% of its total value when explwjui; I « 0.01. Using this as 
a criteriai for infinite length: 


15.131 

where; ■ length of charge within the hot or cold zone for 

■ charge to appear infinitely long. 

If Pe* is small so that eq. [5.5] is satisfied, the characteristic roots 
are given by eq. [5.7]. Substituting cq. [5.7] into eq. [5.13]: 
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Fig. 5.8 ; Effect of cliarg j fXisition within tJie fianace on the axial teinperture 
distri}jution of the charcie. 
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[5.14] 


c, = I 

which is a useful expression for determining . Equation [5.14] is in 
agreement with the experimeital results of Clyne [11] which suggest that 
longer charges and higher Biot numbers tend to stabilize the interface 
position (i.e./ the growth rate is equal to the lowering rate). 

Results from the one-dimensional numerical model of Riquet and 
Durand [9] suggest that the growth and lowering rates are equal when the 
length of charge within ttie hot and cold zones is about 1 . 5 (Bi“l/ 2 ). 
result agrees well with eq. [5.14]; the difference in the constant 
coefficient is attributed ‘ the choice of criteria defining "infinite" 
charge length. 

When the charge lei:^ h is not infinite, appropriate boundary 
conditions must be applied to the eids of the charge. For example, a solid 
pull rod contacting the crucible at the bottom of the charge can be 
approximated by treating it as a simple fin exchanging heat with the 
environment. In this way, a Biot number can be calculated for the end 
boundary condition. Appendix D describes a method to determine the end 
boundary condition for some simple end geometries. The curves of Fig. 5.8 
have been calculated using the same Biot number for the ends of the charge 
as for the circumference. 

CHARGE DIAMETER EFFECT 

Any change in the charge diameter, D, affects the dimensionless 
parameters Bi, ?e, S , , and • Further, changes in D may 

alter the heat transfer coefficients between the furnace and charge as the 
geometry of the furnace cavity changes. The effect of a change in D is 
assessed by reevaluating the necessary parameters and then using the 
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appropriate expressions for the axial gradient. In gaieral, increases in D 
may either increase or decrease the axial temperature gradient in the 
gradient zone depending on the corresponding changes in the other 
parameters. 

As an illustrative example^ consider a symmetric system of infinite 
charge length. The appropriate expression for the axial temperature 
gradient in the gradient zone is given by eq. [5.8]. Since, however, 
changes in diameter are to be analyzed, it is more informative to compare 
the axial gradients of eg. [5.8] based on the dimensional axial 
coordinate, Z: 


= ^ zl >5.151 

Defining the variable: 

^ " iD 


Taking Lq as constant, eg. [5.15] shows that the dimensional axial 
gradient, gL, increases with D when f is negative and vice-versa. 

For an increase in D accompanied by a proportional increase in Dcrf 
(i.e., 8 remains constant), egs. [4.7] and [4.14] show that Bi** increases 
at most linearly with D (assuming that all other thermal parameters are 
constant). Consequently (using eq. [5.16]), f is positive and the 
dimensional axial gradient, gL, according to eg. [5.15], decreases as D 
increases. 

The crucible diameter may, on the other hand, remain constant as D 
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increases. This occurs when increasing the charge diameter within a 
constant diameter furnace cavity at the expense of the thickness of the 
crucible. Fig. 5.9 shows the corresponding behavior of 9 for various Bi 
and K with 5 • 1.25. Accordingly, f is either positive or negative 
depending on the values of Bi and K? the axial gradient of eg. [5.15] may 
therefore increase or decrease with an increase in D. 

In summary, the effect of a change in charge diameter on the axial 
temperature gradi^t at the growth interface should be examined on a case 
by case basis. 
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^ ^ ^ Dcr= constant 



Fig. 5.9 : 'Hie riarameter V (eci,(5.16]) for a syirmetric system, 6 = 1.25 and 

constant 
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Ig£ 7^IMENSIQNftL THERMAL MODEL 

Thermal criteria for satisfactory Bridgman crystal growth include 
requirements for the shape of the melting temperature isotherm (c.f.^ 
section 2.1.2). The one-dimensional model of Chapter 4 does not permit 
evaluation of the isotherm shape since it neglects radial temperature 
variations in the charge. Therefore, a two-dimensional model of the 
Bridgman growth system is developed in this chapter. Results from the two- 
dimensional model are used in Chapter 7 in order to establish the causes of 
radial temperature variations within the charge and to study their 
depeidence on the thermal parameters of the Bridgman growth system. 

THE CCMBTORIC FIN APPROACH 

Hie present two-dimensional thermal model of vertical Bridgman growth 
treats the entire system as a set of coaxial fins. Radial temperature 
variations are therefore accounted for by differences in the fin 
temperatures at the same axial location. The model can thus be considered 
a hybrid between a continuous description in the axial direction and a 
finite difference description in the radial direction. This approach, 
rather than a complete finite difference or finite element formulation, is 
taken for several reasons; 

(1) The primary objective of the two-dimensional model is to 
develop insight into the origin and behavior of radial 
temperature gradients near the growth interface for the 
purpose of controlling the shape of this interface through 
system design. The increased accuracy achievable using a 
fine grid network in a finite difference or finite element 
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model is not considered essential for this objective and, 
furthermore, does not appear warranted for systems for which 
reliable thermal property data are not yet available. The 
present model incorporates all the criti.'al thermal effects; 
it is expected, therefore, to be able to predict the radial 
gradient behavior at least to first order approximations, 
(By utilising enough coaxial fins, the present model should, 
in principle, be as accurate as a finite difference model.) 

(2) The present model is easily implemented on the computer and 
is a simple extension of the one-dimensional moving fin model 
developed in Chapter 4. 

(3) Even with tne simplifications inherent to the concentric fin 
approach, the ^dded complexity of a two-dimensional 
formulatiai necessitates the use of a computer to calculate 
the two-dimensional temperature distribution of the charge. 
The concentric fin model can be simplified, however, without 
loss of pertinent information, to a point which permits 
presentation of some relevant results in analytical form. 


^ I TTE DESCRIBING EQUATIONS FOR THE CD] 




RIC FIN MODEL 


The factors of concern for the development of the concentric fin model 
of the Bridgman growth system depicted in Fig. 2.2 are indicated in 
Figs. 6.1 and 6,2. The two-dimensional thermal model incorporates the 
following assumptions in addition to those described in Section 2,2.2: 


(1) Radial temperature variations within the charge are accounted 
for by modeling the charge as two radially distinct fins. 


Hot Zone 


Gradient Zone 


Cold Zone 


inner charge fin 


solid 


outer charge fin 
crucible 


hot heat pipe 


v////m 


Gradient Zone 
Annulus 


cold heat pipe 


Fiq. 6 . 1 : Schc^^tic of ti\e charge, cwcible, gradient zoikj annulus and 

furnace lieat j)ifX2s in tlie concentric fin model. 
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Ttie inner fin is a cylinder or diameter and has an axial 
temperature distribution The outer fin is an annulus 

of inner and outer diameters* and D* and has an axial 
temperature distribution The diameter is chosen 

so that the inner and outer fins have equal cross sectional 
area. 'Rie temperature of each charge fin is not considered a 
function of radius. The difference in temperature between 
Tjji and Tout indicative of radial temperature variations 
within the charge. 

(2) The crucible is modeled as a coaxial annular fin in contact 
with the outer charge fin. The crucible fin temperature* 
"^cr (Z), is not a function of radius. 

(3) As opposed to the corresponding assumption in the one- 
dimensicxial model* the gradient zone is not adiabatic; i.e.* 
hQ # 0. This assumption is relaxed in the two-dimensional 
model so that the effect of radial heat transfer between a 
non-perfectly insulating gradient zone and the charge can be 
studied. 

(4) The portion of the furnace in the gradient zone is considered 
to be an annulus of conductivity l <<3 and cross sectional area 
A (5 located coaxially with the charge and crucible. The 
temperature of the gradient zone annulus is assumed constant 
in the radial direction; its axial temperature distribution, 
Tq(Z), depends on its thermal interaction with the charge ond 
crucible. Temperature boundary conditions are used at the 
ends of the gradient zone annulus and no heat loss is assum.ed 
at its outer diameter. (This model of the gradient zone 
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annulus was also used in section 4.4.4 in order to determine 
the effect of a non-perfectly insulating gradient zone on the 
axial temperature behavior of the charge.) 

On the basis of the above assumptions, each of the concentric fins is 
analogous to the moving thin tod treated by Catslaw and Jaeger 129)! 
eg. (4.1) (neglecting the transient term by virtue of the quasi-steady 
assumption) therefore describes the axial temperature distribution within 
each fin. In the concentric fin model, the q’ term of eq. t4.ll accounts 
for the heat transfer between radially adjacent fins. In Ajpendlx, P, this 
radial heat transfer Is approximated resulting In the following set of 
coupled ordinary differential equations for the axial temperature 

distribution of the fins: 



[6.1a] 


[6.1b] 


_ p + A 3 + A 4 . ^ 





[6.1c] 


[6. Id] 


where: 
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where: 


A'll ■ P- ) 


«-«’»« Q^Aur? 


Pei.c’ 

Pe.,, = VU(fCf/k)t, 

The coefficients A^, A2r A3, A4, and A5 depend on the radial thermal 
resistance between the concentric fins. They are given by: 




[6.2b] 


[6.2c] 


[6. 2d] 


[6.2e] 


The fin equations for and have been added to and suotracted 
from one another in Appendix F in order to produce eqs. [6.1a] and [6.1b] 
for (p^ and 0^ respectively. Since the cross ' »ctional areas of the inner 
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and outer charge fins are equal, <p^ is interpreted as the mean cross 
sectional temperature of the charge; is a measure of the radial 
temperature difference within the charge. Within the hot and cold zones, 
the furnace temperature, , is the temperature of the hot and cold heat 
pipes, respectively, emd eq. [6. Id] for the gradient zone annulus in not 
needed. Within the gradient zone, ^ is identified with <^(0 
eq. [6. Id] has no Peclet term since the gradient zone annulus is 
staticnary. 


iu SQuniQiiQg ms 




The set of eqs. [6d] has constant coeffici '.r.ts if a separate set of 
equations is ■•mployed for each axial regiOT of uniform thermal properties 
and furnace boundary conditions. As in the one-dimensional model, four 
such axial regions are present in the concentric fin model; the hot and 
cold zcxies and the liquid and solid portions of the gradient zone. (It is 
assumed that the growth interface is in the gradient zone.) The hot and 
cold zones are therel “e described by eqs. [6.1a], [6.1b] and [6.1c]; 
the two regions within the gradient zone include, in addition, eq. [6.1d] , 
representing the gradient zone annulus. 

The solution of eqs. [6.1] within a particular axial region is 
described in many texts (e.g., [36,3"^] U The system of p second order 
differential equations can be changed to an equivalent system uf 2p first 
order differential equations by the definition; 



[6.31 


Performing this transformation and putting the result in matrix notation. 
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PX- P 


[S.41 


where: 


r 1 

<P. 

<0cr 

0cr 

W m 


(hot and cold zones) 


’ ^ 
0c< 



^cr 


(gradient zone) 


(Variables in boldface denote matrices.) The elements of P and F are found 
by direct substitution. In the gradient zone, P is an empty matrix 
The solution to eq. [6.41 is of the form: 


X-C.E 


. + c 




»X, 


[6.5] 


The values of w are commonly called the eigenvalues of the matrix P and the 
values of E are commonly called the eigaivectors of P. The eigenvalues and 
eigenvectors are easily found by standard computer subroutines. The 
coefficients denoted by c in eq. [6.51 are the unknown constants of 
integration. Xp denotes the particular solution of eq. [6.4]. 

Equation [6,5] is written for each axial region. The hot and cold 
zones each contribute six unknown constants of integration; the liquid and 
solid parts of the gradient zone each contribute eight unknown constants of 
integration, ^ncre are, therefore, 28 unknown constants of integration 
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which must be determined by using boundary conditions of (a) equality of 
temperature and (b) continuity of flux between corresponding fins of 
adjacent axial regions as well as (c) heat transfer coefficients at the 
ends of charges of finite length and (d) imposed temperature boundary 
conditions at the ends of the gradient zone amnulus. Application of the 
boundary conditions results in a set of 28 simultaneous linear algebraic 
equations whose unknowns are the required constants of integration. The 
solution of these simultaneous equatiois is easily accomplished by standard 
computer subroutines. 


JLI APPRQXIMAnC3N FOR THE 




RFACE BQCNDARY OONDmON 


The solution of the fin equations requires a trial and error procedure 
when the extent of each axial region is not initially known. For example, 
the one-dimensional model of Chapter 4 shows that the interface location, 
, and the interface temperature, 9i , are functionally related. When 
9i is chor-»n as the indepaident variable, must be found by iteration 
(c.f., the analytical results of Appendix B). When is taken as the 
independent variable, the location of the interface boundary conditions is 
known a priori resulting in a direct solution for 9^ as well as the entire 
axial temperature distributicxi. 

Figure 6.3 shows that the axial position of the growth interface in 
the inner and outer charge fins is, in general, not the same; they are the 
same only when the interface is flat. The correct application of the 
concentric fin model requires that a separate set of eqs. [6.1] be written 
for the short region between the axial locations ot the interface. In this 
case, the interface temperature, 9i , depends on the interface locations of 
both the inner and outer charge fins, C,- ■ and Q- ^ respectively: 

Lj iff L ^ OU w 
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Approximated Interface 
Location 



Location 


Fia. 6.3 ; ADoroxiination for the interface boundar 
conditions used in the concentric fin model. 
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f ( 9 - r. . r. ) ^ 0 t6.6] 

Since this functional relationship is not known a priori, it is not 
possible to choose initial values of and ^c^out which satisfy 

eq. [6.6]; the solution therefore requires an iterative procedure. The 
present section describes an approximate model of the interface which 
precludes the use of a trial and error solution for the concentric fin 
model. 

The length of the region separating the interface in the inner and 
outer charge fins is normally small since it is desired to -have an 
interface of small curvature. The preset model neglects this axial region 
and applies the boundary conditions of the interface at a single axial 
location for each of the inner and outer charge fins. The axial location 
of the interface is taken to be the location where the mean charge 
temperature, ^ , is equal to the interface temperature, 9^ (see 
Fig. 6.3). In this sense, the axial location where ^ 9^ denotes an 
"average" interface position (i.e., between the actual interface positions 
of the inner and outer charge fins). The value of at this location, 
, is indicative of the curvature of the interface. 

The use of the above approximation for the interface boundary 
condition permits the a priori choice of the "average' interface location 
and, therefore, the solution for the temperature field in the charge is not 
a trial and error procedure. The approximation introduces no error if 
Rr = 1 and PesRu = 0 or if the interface is flat. Otherwise, when the 
interface curvature is small, the error should likewise be small and the 
predictions or the present model concerning the radial temperature 
variations in the charge are expected to be correct to a first-order 
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approximation. 

The interface flux boundary conditions, assuming that the melt/crystal 
conductivity change and the generation of latent heat occur at the same 
axial location in each of the charge fins, are; 


AC, '$ 


[6.7a] 



[6.7b] 


[6.7c] 


[6.7d] 


where the subscripts "L" and "S" refer to the melt and crystal side of the 
interface, respectively. Equations [6.7] assume that the growth rate, R, 
is equal to the lowering rate, V (see section 4.4.2). 




RS QF THE QON 




FIM MODEL 


The thermal parameters required by the concentric fin model are those 
necessary for the one-dimensional moving fin model (see Table 4.1) and, 
additionally, those listed in Table 6.1. 


CDMEARISQN WITH THE TTO-DIMENSIQNAL MODEL OF FU AND WILCDX [261 
The figures of Fu and Wilcox [26] plot the shape of various isotherms 
within the charge by the use of two-dimensional finite difference computer 
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parameter 

definition 

cemnents 

BiG,s 

h^D/kg 

BiG,L • 




Peer 

VD( f Cp/k)cr 

Peer • Pes/Ro< ,s 

\ i-xjz) 







Table 6.1: 'Ihemal parameters of the concentric fin model, in addition 
to those listed in Table 4.1. 


solutions. Isotherm shapes taken from their paper are used in this section 
in order to check the results of the present two-<3im«isional model. 

In the figures which follow, the shape of an isotherm is indicated by 
the axial distance between the location of ths isotherm at the center of 
the charge and the outer diameter of the charge and is denoted by ^ ^ . 

is scaled directly from the isotherm plots of Fu and Wilcox [261; the 
representative axial location of each isotherm is taken to be the axial 
location of the isotherm at the diameter (c.f.. Fig. 6.2). The preseit 
model predicts and as functions of C, . In Appendix G, r,(p is 
related to and through eqs. [G.2] and [G.5] : 

r = ^ (6.81 

The system modeled by Fu and Wilcox [26] has an adiabatic gradient 
zone, infinite charge length, and does not include a crucible, generation 
of latent heat, or change in thermal conductivity at the growth interface. 
The results of the concentric fin ir.odel are therefore obtained using the 
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following parameter values: ^ ■ 1> R{^ > Ir Rh * Or and an adiabatic 

gradient zone. Four cases, each with different Biot or Peclet numbers, are 
compared in the graphs of Fig. 6.4; the parameters tested are listed in 
Talkie 6.2. It is seen that there is good agreement between the results of 
the two models. (The apparent scatter in the results taken from Fu and 
Wilcox [26] is attributed to difficulty in scaling the values of from 
their plotted isotherms.) It is therefore concluded that the present 
concentric fin model of the charge predicts, at least to first-order 
approximation, the correct two-dim^sional thermal behavior of the charge. 

test Of m EFFECTIYE BIOT MIBER 

The effective Biot number, Bi**, is developed in Sections 4.2.2 and 
4.2.3 in order to account for the effects of crucible conductivity and 
thickness as well as radial temperature variations within the charge on the 
axial temperature distribution of the charge. Section 4.6 provides a 
partial verification of the effective Biot number concept by comparing the 
one-dimensional results of this work, using Bi**, to the results of the 
two-dimensional model of Fu and Wilcox [261. This comparison, however, 
included only the effects of radial temperature gradients within the charge 
on Bi** since the model of Fu and Wilcox [261 does not provide for a 
crucible. The present two-dimensional model can be used to verify the 
quantitative accuracy of Bi**, including the effects of the crucible. For 
this purpose, a symmetric system with Bi <3 ■ 0 and of infinite length is 
used. The applicable relation for the axial temperature gradient in the 
gradient zone from the one-dimensional results of C3Tapter 5 is eq. [5.81: 



original 
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Adiabatic gradient 

zone: 

BiQ«10'4 

Rs«IO-« 

1 

Symmetric system 

1 

Xg»I.O 

A ‘ 


\ i 

concentric fin model 

\ 

0 X Fu and Wilcox 

\ • 


hot zone- 


gradient zone 


Fig. 6.4a : Corparison of isotherm shapes frcn the results of the 

concentric fin model and Fu and Wilcox [26]; Bi = 0.4, 4.0, and 
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Fia. 6.4b: Cbrparison of isotherm sha^s fron the results of the^ 
concentric fin model and Fu and Wilcox [26] : Bij^ = 2.0, Bi(; - 4.0 
and Pe = 0.0. 
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i-ig. 6.4c; Comparison of isotherm shares from the results of the 
concentric fin model and Fu and Wilco:: [26] : Bi = 2 0 and 

Pe = 0.4. 
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Figure from 
Fu & Wilcox [26] 

— 

Bifl 

■ 

Bic 



Pe 

Fig. 6.4a 

Fig. 4 

0.4 

0.4 

0.0 

Fig. 4 

— !■ 

4.0 

0.0 

Fig. 6.4b 

Fig. 5 

BI 

4.0 

0.0 

Fig. 6.4c 

Fig. 7 

m 

■a 

0.4 


- 1.0 S - 1.0 

■ Ac ■ % ■ 0.0 

R|^ ■ 1.0 

adiabatic gradient zone 

Table 6.2: Parameters used for comparing the results of 

Fu and Wilcox [261 and the present concentric 
fin model. 


The system used in this comparison is further specified by taking » 1 
and <S » 1.25. 

The axial gradient given by eq. 15.81 is compared in Table 6.3 to that 
determined by the concentric fin model for several values of Bi and K. It 
is seen that the agreement is excellent. It is concluded, therefore, that 
Bi** is an accurate approximation for the effects of the crucible on the 
heat exchange between the furnace and the charge. 


siGENvamss or the 




MC FIN ?CPa 


The eigenvalues reflect, in a qualitative sense, the thermal 
interactions among the fins which cause their temperatures to equilibrate 
in the axial direction. Larger eigenvalues suggest faster equilibration. 
It is expected that the eigenvalues increase as the ratio of axial to 
radial thermal resistance increases. As an illustrative example, consider 
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orie^imens ional 
mouel 

two-dimensional 

model 

Bi 

K 

Bi** 

eqs. [4.7] 
and [4.14] 

dVdd; 
eq. *5.8] 




0.473 

-0.407 

-0.415 


D 

0.684 

-0.453 

-0.459 

1,0 


0.660 

-0.448 

-0.453 


3.0 

0.427 

-0.395 

-0.397 


10.0 

0.183 

-0.300 

-0.299 


0,1 

0.103 

-0.243 

-0.244 

1 


0.3 

0.101 

-0.241 

-0.242 

0.1 

1.0 

0.078 

-0.219 

-0.219 


3.0 

0.046 

-0.177 

-0.176 


10.0 

0.019 

-0.120 

-0.120 


S • 1.25 ?e « 0.0 
% - 1.0 » 1.0 
% » 0.0 

adiabatic gradient zone 


Table Testing the effect of the crucible on 

Bi** by using the concentric fin model. 
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the r.ovinq fin model of Chapter 4: tne Biot munber expresses the ratio of 
axial to radial thermal resistance in this case. Equation [5.8] and 
Fig. 5.2 show that as the Biot number increases, the significant 
temperature difference in the model, , decreases more rapidly with 
distance into the hot and cold zones. 

As eq. (6.51 indicates, the eigenvalues only partially express the 
thermal behavior of the concentric fins. The eigenvectors, B, and the 
constants of integration, c, combine to determine the importance which ©ich 
eigenvalue has for each of the fin temperatures. Nevertheless, a study of 
the parametric behavior of the eigenvalues serves to emphasize the nature 
of the thermal interactions occur ing among the concentric fins. 


6.8.1 Hot and Cold Zones 

By substituting a solution of the form: 


\<Pr.' 


1 

i 

U 


• 

Ca 1 

[<f>cr 


1 


[6,9] 


into eqs. [6.1], the characteristic equation for the eigenvalues of the hot 
and cold zones is found to be; 
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The six roots cf equation (6,10] are the eigenvalues of the hot and cold 
zones. By inspection of eqs. (6.21 which define the coefficients A2» 
A3, A4 and A5, it is seen that the eigenvalues depend on four parameters: 

S * ^loc» ®^loc ^®loc* 

Wlien P«ioc ■ Of ((5.10] becomes a cubic polynomial in w2 indicating 
that the eigeivalues, in this case, appear in pairs of opposite sign. When 
Peioc ^ 0 but K]_oc “ If - t may be shown by direct substitution that the 
eigenvalue pairs have the following form: 



(6.11a] 

I 

II 
1 

(6.11b] 

vdiere: a is a positive number 


Whei Kxoc ^ ^f results suggest that the eigenvalue pairs are 
number whose magnitude is of the order Pexoc^Kioc* 

separated by a 


(6.12a] 

w - -a 

(6.12b] 

where: t. is of order unity 


Since the Peclet number is typically small, it is appropriate to study the 
behavior of the eigenvalue pairs assuming that Pe » 0? in this way, the 
order of the characteristic equation, eq, (6.10], is halved. 

Figure 6.5 illustrates the variation of the three eigenvalue pairs as 
a function of S with and Bi^Qc parameters. The lowest set of 
curves is seen to vary greatly with Bi^oc* Further, the values of the 


eigenvalues represented by the lowest set of curves is found to be closely 


- 124 - 


Eigenvalues for the hot and cold zones, Pe 







ORIGINAL PAGE IS 
OF POOR QUALITY 



0.10 1.0 


(«-l) 

. 6.5b ; Eigenvalues for the hot and cold zones, Pe = 0.0 
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predicted by eq. (5.7] : 


original 

OF PoQff 

(5.7] 

(The largest difference between the eigenvalues predicted by eq. (5.7] and 
those in Fig. 6.5 is 261 at S ■ 5 and Biio^ * 

difference decreases for smaller S ai^ Bi.) These eigenvalues, therefore, 
can be associated with the variation of the mean charge temperature and are 
denoted by w^. As w^ increases, tends to approach more rapidly the 
temperature bomdary conditions of the furnace. 

The eigenvalues represented by the middle and upper set of curves in 
Fig. 6.5, d«ioted by w^ and W 2 respectively, are relatively insensitive to 
Biioc« They do exhibit, however, a strong sensitivity to and S , 
parameters which describe the thermal interaction between the charge and 
the crucible. 

For the limiting case of i 1, i.e., no crucible, it is seen that 
wj approaches a constant value (approximately 7.6) independent of S 
Kioc* This value is therefore most easily identified with the thermal 
characteristics of the charge alone. Results derived in section 7.2 
indicate that radial temperature gradients in a charge without a crucible 
attenuate at a rate proportional to exp(-7.63C ). Thus, w^, is indicative 
of the rate of decay of for systems when & is small. The value of W 2 
is very much larger than either wj^ or w^ for small values of & suggesting 
that the crucible temperature equilibrates rapidly with the temperature of 
the otter charge fin. 

A second limiting condition is that of large S . Figure 6.5 shows 
that W 2 approaches the value of 7.6 as S increases; further, the approach 


Quality 
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to the limiting value takes place at smaller values of 6 as Kiqc 
decreases. These facts suggest that the charge is, in this case, 
"insulated" in the sense that heat transfer within the charge is much 
better than that between the charge and the crucible. 


6.8.2 

The characteristic equatiai for the eigenvalues in the gradient zone 
is determined in a manner analogous to that used to develop eq. [6.10]: 


Kioc ' 




As) ~ ^ 


[6.13] 


L ^\oc. 

Equation [6.13] is a polynomial of eighth order in w; there are therefore 

eight eigenvalues in the gradient zone. 

In eq. [6.13], the term of lowest order in w is the linear term: 


-W 



- 


[6.14] 


Since there is no term of order zero in w, one of the eigenvalues is zero. 
Further, when Peioc = shows that the linear term disapijears; 
in this case, a second eigenvalue is zero. The solution for the fin 
temperatures given by eq. [6.5] is valid only when there are no repeated 
roots. When Peioc * Of the term which would include the second zero 
eigenvalue is replaced by a term linear in C . 
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The linear component of equation [6.5] in the gradient zone can be 
demonstrated by multiplying eq. [6,1a] by eq. [6.1c] by kc^Acr and 
eq. [6.1d] by 1C(3 Ag and then adding. With Pe^oc ' the result is: 





= 0 


[6.15] 



[6.16] 


0 « 

^lot A Ac/> Afc 

^ is interpreted as an average temperature of the 
charge/crucible/gradiait zone annulus combination. Equation [6.15] results 
in a linear variation in <p within the gradient zone. The physical 
interpretation is that there is only internal heat exchange between the 
fins of the gradient zone; there is no external heat addition or removal to 
the charge/crucible/gradient zone annulus combination by virtue of the 
assumption that the gradient zone annulus is adiabatic at its outer 
circumference. Therefore, neglecting Pe, the axial conduction of heat is 
constant within the gradient zone. The definition of gi as an average of 
t and weighted according to their respective axial heat 
conducting capacities is demonstrated by eq. [6,15] to correctly express 
the constant total axial heat transfer. 

As occurs for the eigenvadues of the hot and cold zones, results show 
that the eigenvalues of the gradient zone occur in pairs of opposite sign 
which are separated by a number whose order of magnitude is that of 
Peioc^^loc (c.f., eq. [6.12]). The partner of the zero eigenvalue 
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therefore has a value of order of magnitude Peiog. Since Pe^oc is 
typically small, the exponential function of eq. [6.5] with this eigenvalue 
is very nearly linear in the gradient zone. 

The additional fourth pair of eigenvalues introduced by the gradient 
zone annulus within the gradient zone are associated with the rate of 
equilibration of ^ with ^ . An analytical expression for this 
eigenvalue based on a symmetric, one-dimensional system was developed in 
section 4.4.4: 



Table 6.4 compares the values of wq determined from the concentric fin 
model and from eq. [6.181. The comparison is favorable, indicating that 
this eigenvalue pair of the conceitric fin model is correctly associated 
with the rate of equilibration of dj. and . 

The eigenvalue pairs w^ and W 2 are shown in Section 6.8.1 to be 
relatively insensitive to conditions external to the charge and crucible. 
Table 6.5 demonstrates that the additional parameter introduced into 
eqs. [6.1] by the gradient zone annulus, i.e., Rq, also has a second order 
effect on these eigenvalues. Therefore, approximations for w^ and W 2 of 
the gradient zone may be obtained from Fig. 6.5. 
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*ks,loc 



0.01 



10.0 


0.01 


.7908 

.6633 

.2943 

.2828 

.1810 

.2098 

Bi<3 

0.1 

7.729 

6.356 

2.476 

2.098 

m 

.5663 

.6633 

1.0 

23.22 

21.10 

7.601 

6.63 



2.653 

2.828 

1.644 

2.098 


10.0 

43.34 

44.94 

m 

3.967 

6.325 

2.742 

4.690 



^ from COTcentric fin model 
WQ fran eg. [6.18] 


Pes ■ 0.001 Kiqc ■ 1.0 % ■ 1.0 

Table 6.4; Comparison of wq determined from the concentric 
fin model and from eq. [6.18] 
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Bi<3,loc 


I^le 6.5: 





Bg ,1 oc 

0.01 

0.1 


10.0 

0.01 

9.635 

4.838 

9.635 

4.837 

9.635 

4.837 

9.635 

4.837 

0.1 

9.641 

4.824 

9.637 

4.859 

9.637 

4.854 

9.637 

4.853 

1.0 

9.631 

4.828 

9.686 

4.729 

9.657 

5.054 

9.656 

5.009 

10.0 

9.631 

4.827 

9.577 

4.763 

9.741 

6.255 

9.718 

5.569 


key: 


wi 

W2 


Pes » 0.001 Kioc * 1.0 % * 1.0 

The relationship of and W2 in the gradient zone 
to the parameters BIq^ioc ^nd Rg,1oc* 
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2^ RADIAL TEMPER ATURE VARIATIONS WITHIN THE CHARGE 

Radial temperature variations near the growth interface are shown in 
this chapter to result from several effects usually present in Bridgman 
crystal growth. Systems with an adiabatic gradient zone, first studied by 
Chang and Wilcox [17]/ addresses only one of these factors. The other 
factors have not been adequately considered in the literature: non- 
perfectly insulating gradient, zone (he ^ 0), and thermal conductivity 
change and/or generation of latent heat at the growth interface in the 
presence of a crucible. Computer solutions of the concentric fin model are 
used to graphically illustrate the behavior of the radial temperature 
variations under the influence of these factors. Special attention is 
given, however, to the developmeit of approximate analytical expressions 
which serve to emphasize the functional dependence of the radial 
temperature variations on the various system parameters. 

The variable defined by the concentric fin model of Chapter 6, is 
used to indicate the radial temperature variations within the charge. 
The results for to be shown in the following sections are convaiiently 
approximated to the interface shape by the following relation developed in 
Appendix G: 

where M is the radius of curvature of an isotherm in number of charge radii 
at a particular axial location. (N is positive for an interface which is 
concave toward the crystal as shown in Fig. 2.1.) Accordingly, the degree 
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of flatness of isotherms is enhanced by large axial temperature gradients 
and sma^l radial temperature gradients 


2A GftUSES QF RftDIAL TEMPEKflURE VaBIAIIQMS WIIHIM TBE GRftCian: 2QME 

Subtracting eq. [6.1b] from eq. [6.1a ] , and assuming that the Peclet 
term is negligible (c.f., section 5.1) , yields: 







[7.2] 


The particular solution of eq. [7.2], denoted by , may be solved by using 
the appropriate Green's function: 


exp(-A';‘U-xO 


Using eq. [7.3], the particvdar solution of eq. [7.2] is; 



\ 

\S-.2t» 





[7.4] 


where the value of A2 given by eq. [6.2b] has been used, x is a dummy 
variable of integration, and and Ch are the axial locations of the 
cold and hot ends of the charge, respectively. (Nc ce that eq. [7.4] is not 
the complete solution for since and are not independent by 
virtue of the coupled nature of the fin equations. Nevertheless, the 
actual distributions for and <p^ must satisfy eq. [7.4].) 

is non-zero when the combined axial heat conduction within 
the inner and outer charge fins is not constant with axial position; i.e., 
when there is heat addition to or removal from the charge. Equation [7.4] , 
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therefore, clearly demonstrates that radial temperature gradients within 
the charge are caused by heat exchange between the charge and its 
surroundings. Such heat exchange occurs: 


1) if, in the idealized case of no crucible, there is a 
differmce between the temperatures of the outer charge fin 
and the furnace, and the heat tr^’^sfer coefficient between 
them is not zero. 

(2) if, when a crucible is pres«it, there is a difference between 
the temperatures of the crucible and the outer charge fin. 

Equation [7.4] expresses <p^ at a particular axial location as the 
superposed effects of heat transfer to the charge along its aitire length. 
However, the effect of heat transfer to the charge at one location, x, on 
the radial temperature gradient at another location, C , attenuates rapidly 
with the distance between the two locations, being proportional to the 
factor exp(-7.63l C - xl). For example, when I C; - xl « 1, the value of 
this exponential function is 0.00049. "nierefore, the radial variation in 
temperature at the location C is primarily affected by heat exchange with 
the charge within the range ^ f 1 (i.e., one charge diameter to either 

side of the giv«^ location). 


7.1.1 Analysis of Radial Tgnperatu r p Variations In the Gradient Zcaifi 

The total solution for <P^(C) includes, in addition to eq. 17.4), the 
homogeneous part of the solution; 

r+s fr\ r r ^ 17.51 
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Equation [7.5] is applied separately to each region of the charge of 
uniform properties; given the assumptions of the concentric fin models 
there are two such regions — the liquid and the solid portions of the 
charge. The coefficients and C2 in eq, [7.51 are found by applying 
boundary conditions at the ends of the charge and at the growth interface. 

Radial temperature variations near the growth interface are of primary 
concern. Since the interface is assumed to be in the gradient zoner 
eq. [7.5] is applied, in the remainder of this chapter, only to the 
gradient zone; i.e., X^/2 . It is also assumed that there is 
at least one dieuneter of chirge in both the hot and cold zones; 

Ch ‘ l') !7.6al 

Cc C7.6b) 

As noted earlier, heat exchange with the charge at distances greater than 
about one diameter from the ends of the gradient zone contribute little to 
the radial temperature variations in the gradient zone. Equations [7.61 
thus permit the limits of integration in eq. [7.5] to be contracted if 
desired, and fuiuher, the effects of the boundary conditiats at the ends of 
the charge on can be neglected. 

Wlien the solid and liquid thermal properties of the charge are equal, 
i.e., Rr ■ 1, eq. [7.5] need be applied only once to the entire length of 
the charge. The homogeneous part of the solution disappears under the 
restrictions of eq. [7.6] and the total solution for within the 
gradienc zone becomes; 
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where: r ^ r ^ k i 

Cl s C s XJz 


The variables Il(? ) and Is(C ) represent the respective cor*tributions of 
radial heat exchange with the liquid and solid portions of the charge to 
the value of at any axial location ^ within the gradient zone. They 
are defined as follows: 
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[7.11a] 


[7.11b] 


The first term in each of eqs, [7.9] and [7,10] represents a 
contribution to which is produced at the growth interface and which 
attenuates with distance from the interface. This contributiai is absent 
when % = 1. 

The value of at the interface, , is determined from either 
eq, [7,9] or eq, [7,10]: 


Kk"*" I 

An attractive possibility for the quantitative utilization of eqs, 
[7,9] and [7,10] is to provide from the results of one- 
dimensional modeling; that is, to approximate with . In 

this way, a two-dimensional temperature distribution is obtained only 
through these equations and a one-dimensional model. The error in 
calculated in this manner would be due to approximations in (1) the one- 
dimensional model which estimates and (2) the concentric fin 

model which provides eq, [7,2], 

One-dimens icnal models do provide good approximations for the charge 
mean temperature distribtuion, • Taking two derivatives of such a 

distribution in order to obtain a good approximation for is more 

likely to be unreliable. Results presented later in this chapter, in fact. 
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demonstrate that the presence of a crucible produces a distribution of 
radial heat exchange with the charge which is not predicted by a one- 
dimensional model; therefore/ the use of calculated from one- 
dimensional models is not appropriate for use in eqs. [7.11]. 
Nevertheless/ ev«» in these caseS/ the Green's function approach developed 
in this section may be used to obtain a qualitative understanding of the 
behavior of radial temperature variatiais within the charge. 

PENEmftEQM Qf RADIAL TagEBflmRE VaRimONS. TBE G RflDia n: ZQtffi 

FROM THE HPT AtP QDLD ZONES 

Heat transfer to the charge in the hot zone and from the charge in the 
cold zone must occur so that the required axial temperature gradient at the 
growth interface is obtained. Such heat exchange, according to the results 
of section 7.1, produces radial temperature variatiais within the charge 
which, even within an adiabatic gradient zone, influoice the shape of the 
growth interface. 

Fu and Wilcox [26] studied the penetration of radial temperature 
variations into an adiabatic gradient zone for systems without a crucible 
through the use of a two-dimensional finite difference model. Isotherm 
shapes indicate, for various Bi and Pe, that radial temperature variations 
within the charge in the hot and cold zones attenuate rapidly from the ends 
of the gradient zone. If the gradient zone is sufficiently long, there 
exists a region within which the shape of isotherms is quite flat. 

1.2j 1 Systems Without a Crucible 

Figure 7.1 shows the variation of in an adiabatic gradient zone, 
obtained by the concentric fin model of Chapter 6 for systems without a 
crucible. (Note that the vertical sca..e of Fig. 7.1 changes from 
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hot end 


^ in gradient zone 


cold end 


Fig . 7.1 : Variation of in an adiabatic gradient 

zone vdthout a crucible. The system is s'_mnietric . 
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logarithmic to linaar at 1 so that both nagativt and positiv* valuta 
of may be plotted in the same figure.) Accordingly, it is seen that 
attenuates at a rate approximately proportional to exp(-7.6 ) where 

s* is measured from the ends of the gradient zone. It is also seen that 
in the gradient zone increases with the Biot number of the hot and cold 
zones reflecting greater amounts of radial heat transfer to the charge in 
these zones. 

The variation of the radius of curvature of isotherms for the same 
cases considered in Fig. 7.1 is shown in Fig. 7.2. These curves were 
calculated using eq. [7.1] along with values of 0^ and determined 

from the concentric fin model. Accordingly, there is a central regic»i of 
the gradient zone where the isotherm shapes are quite flat. To obtain an 
interface which is slightly concave toward the crystal, i.e., small but 
positive N'i, requires that the interface be located slightly to the hot 
side of the cent*?r of the gradient zone. Also seen in Fig. 7.2 is that 
iscjtherms are more curved as Bi increases. Therefore, for the parameter 
values considered, increasing the Biot number of the hot ind cold zones has 
a greater effect on increasing than upon increasing the axial 
temperature gradient of the charge in the gradient zone. 

The effect of unequal hot and cold zone Biot numbers in an otherwise 
symmetric system is shown in Fig. 7.3. As the hot zone Biot number 
increases relative to the constant cold zone Biot number, becomes more 
positive in the hot portion of the gradient zone. In this region, 
therefore, isotherm shapes become more concave toward the crystal .and the 
location of the flat isotherm shifts slightly toward the cold zone. The 
effect of changing thermal conditions in t)ie hot zone attenuates rapidly 
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Fig. 7.2 ; Variation of isotherm shapes for the same cases 
considered in Figure 7.1. 
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hot end C, in the gradient cold end 

zone 


Fig. 7.3 : The effect of unecual hot and cold zone Biot 

numbers on in an adiabatic gradient zone. The 
s\'’steins are otherwise s^mnetric. 
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from the end of the hot zone and affects only slightly the radial 
temperature variations near the cold end of the gradient zcaie. 

Figure 7.4 shows the variation of in systems which are symmetric 
except Rr # 1 . In the systems considered, the cold zone Biot number 
remains constant; the hot zone Biot number therefore changes with Rr 
(B ig ■ Bic/RR). As Rr increases, it is seen that <p^ in the hot part of 
the gradient zone decreases whereas in the cold part of the gradient 
zone is little affected. This behavior is similar to that of a changing 
value of Big exhibited in Fig. 7.3. Further, a disturbance appears at the 
interface which increases as increases and which dies out with distance 
from the interface. This disturbance represents the first term in eqs. 
[7.9] and [7.10] where it is shown that it appears only when R^ ^ 1. 

An analytical expression for the variation of <p^ in an adiabatic 
gradient zone may be obtained by using determined by the moving fin 
model in order to approximate [7.11]. The charge is 

assumed infinite in leigth; in this case, 8m is described by a simple 
exponential function in both the hot and cold zones. Further, the Peclet 
number is assumed small enough to satisfy eq. [5.6]; the argument for the 
exponential function is then given by Wj,^ in eq. [5.7]. Hence; 

hot zone: ^ ® “ ^g, 

• i ■* 0 ■ [7.13a] 

cold zone: \^ji 

Substituting the second derivatives of eqs. [7.13] into eqs. [7.11] and 
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noting that 


0 in the adiatotic gradient zone, yields: 




5.S2lfr(»i + ZCBl«)''»] 


expl' ?. Cc* X 4 /2)| 


[7.14a] 


Where: - X<a/Z < ^ ^ Xo./2- 

Bquaticn [7.14a] represents the effect of radial heat transfer between the 
charge and the furnace in the hot zone on the radial temperature variations 
in the gradient zone. The effect of the hot zone dies out exponentially 
with distance from the hot end of the gradient zone at approximately the 
same rate observed in Fig. 7.1. Equation [7.14b] has an analogous 
interpretation with regard to heat transfer in the cold zone. 

Substituting eqs. [7.14] into eqs. [7.9] and [7.10] produces the 
desired approximation for in the gradient zone due to attenuation of 

radial gradients from the furnace zones. Evaluating eqs. [7.14] at ^ » 2^,- 
and substituting the result into eq. [7.12] yields at the growth 
interface. 

The one-dimensional results required by eqs. [7.14] are the effective 
Biot numbers and the mean charge temperatures at the hot and cold ends of 
the gradient zone. For a system with Big » 0 and which has a negligible 
Peclet number, the axial temperature gradient in the gradient zone, G, is 
constant except for a step change at the growth interface. The 
temperatures at the ends of the gradient zone are therefore approximated by 
the following expressions: 
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iSi.- 

9m(X<i./Z)= &v«C5;) <^s‘ 


where: Pe satisfies eq. [5.61 

Bi<3 ■ 0 

->a/2 SS; « X<i/2. 


[7.15b] 


The non-dimensional interface temperature, , is given by eq. [B.10], 

Gl is given by eq. [B.91 , and Gs is related to Gl in eq. [4.251. 

The circled points in Figs. 7.1, 7.3, and 7.4 indicate values of 
calculated through the use of eqs. [7.141. The values of Gl, Gs and 9-, 
required by eqs. [7.141 have been calculated us'lr.g eqs. [7.151 as described 
above. It is seen that the values of obtained by the two raethods agree 
very closely indicating that the moving fin model provides a good 
approximation for the systems considered, i.e., those without a 

crucible and with an adiabatic gr-dient zone. The largest difference occurs 
for the largest value of Bi reflecting the slight decrease in accuracy of 
the effective Biot nuiiber as Bi increases. 

7.2.2 ^ Crucible 

The effects of the presence of a crucible on radial temperature 
variations in an adiabatic gradient zone are shown in Figs. 7.5 and 7.6. 
The systems examined are symmetric and infinite in length; Fig. 7.5 
considers variable S (with K - 1) while Fig. 7.6 considers variable K 
(with & » .1.5). The primary effects of the crucible are seen to be a 
reduction in the value of at the ends of the gradient zone and a 
reduction in the rate at which (p^ attenuates toward the center of the 
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Fig. 7,5 : Variation of ^>a in an adiabatic gradient zone 

with a crucible. The systems are syirmetric with 
variable diameter ratio 6 (K = 1.0). 
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hot end C, in gradient zone cold zone 


Fic. 7.6 ; Variation of in an adiabatic gradient zone 
with a crucible. The systems are .'=_'r..Tetr-! c 'vith 
variable thermal conductivity' K (5 = 1.5). 
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gradient zone. 

The reouction in the value of <p^ at the ends of the gradient zone 
wh«i a crucible is present is avcributatU: to the effect of the crucible on 
the effective Biot number, Bi**. In most cases, the presence of a crucible 
decreases Bi** compared to a charge without a crucible (c.f., 
section 4.2.2). A reduction in the heat exchange between the charge and 
furnace is therefore expected, resulting in decreased radial temperature 
variaticxis within the charge, llie individual curves of rigs. 7.5 and 7.6, 
labeled with their respective values of Bi**, demons rate the correlation 
between and Bi**. 

The reduction in the rate of attenuation of cp^ when a crucible is 
present results from radial heat exchange between the crucible and the 
charge near vhe ends of the gradient zone. Such heat transfer occurs when 
the crucible and the charge temperatures are not equal, value of zero 
for the local Biot number is no longer sufficient to eliminate heat 
transfer to the charge in the gradient zone. 

Ibe variation of gradient zone, is 
shewn in Fig. 7.7 for the case of S ■ 2 from Fig. 7.5. It is seen that 
the difference in the crucible and charge temperatures created by heat 
transfer in the hot zone furnace extends into the adiabatic gradient zone. 
The resultant effect cm radial heat exchange with the charge is indicated 
by the distribution of in Fig. 7.8 (lower solid curve). It is 
seen that is not zero in Jie adiabatic gradient zoie; the effect 
of the crucible is to extend radial heat transfer taking place in the hot 
zone into the gradient zone. Analysis of eqs. [7.11] verifies that this 
extension of radial heat exchange serves to decrease the apparent rate of 
attenuation of at the ends of the gradient zone. 
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hot zone . 


.gradient zone 


:Xg/2 


Fig. 7.8: Distributions of and near the hot end of the 

gradient zone, calculated by the moving fin and concentric fin models. 
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quality 

Assuming that the Pe term in eq. [4.17] is negligible, the one- 
dimensional moving fin model of Chapter 4 predicts: 



Since Bi^ » 0, the moving fin model predicts that ■ 0 in an 
adiabatic gradient zone, whether a crucible is present or not. Ihe dashed 
curve of Fig. 7.8 is the distribution of determined from the 
moving fin model. The moving fin model is unable to account for the 
particular nature of radial heat transfer betwe^ bie charge and crucible 
caused by the abrupt change in furnace boundary conditions at the ends of 
the gradient zone. Although the error is of minor consequence in 
calculating the distribution of Qm, results from the moving fin model, in 
particular inadaquate. to describe the local variation of 
radial heat transfer near the ends of an adiabatic gradient zone when a 
crucible is present. 

Table 7.1 compares values of at the end of the hot zone (i.e., at 
^ -0.5) calculated by eq. [7.14a] (i.e„ with the use of the moving 
fin model) and by the caicentric fin model for the cases shown in Figs. 7.5 
and 7.6. It appears that eq. [7.14a] predicts well the qualitative trends 
as K and S change as well as being, to a certain degree, quantitatively 
accurate. This conclusion provides further evidence for the correlation 
between Bi** and noticed in Figs. 7.5 and 7.6. Such a 
favorable comparison may seem, however, surprising in v:'ew of the 
limitations in the moving fin model discussed in the proceeding paragraph. 
This apparent anomoly is explained by a closer examination of the 
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moving fin 

model 

(-Xc/Z) 

Kloc 

S 

Bi** 



eq. [7.14a] 

concentric 
fin model 

1.0 


0.471 

-0.407 

0.703 

0.00407 

0.00407 

1.0 

IH 

0.291 

-0.350 

0.675 

0.00284 

0.00263 

1.0 

2.0 

0.205 

-0.312 

C.656 

0.00216 

0.00193 

1.0 

5.0 

0.0619 

-0.199 

0.600 

0.000797 

0.000683 

0.25 

1.5 

0.357 

-0.374 

0.687 

0.00331 

0.00318 

0,50 

1.5 

0.356 

-0.374 

0.687 

0.00331 

0.00310 

2.0 

1.5 

0.200 

-0.309 

0.654 

0.00212 

0.00197 

4.0 

1.5 

0.121 

-0.258 

0.629 

0.00141 

0.00132 


Bi “ Bic « Bi{j =* 0.5 
Symmetric system 

- 1.0 

Xy * “ ^CD 

liable 7.1: Comparison of eq. [7.141 and the concentric fin model to 

predict <t>a at the end of the gradient zone. The systems are 
those examined in Figs. 7.5 and 7.6 


distributions of and shown in Fig. 7.8. 

The abrupt change in predicted by the moving fin model at the 

end of the gradient zone becomes smoothed out in the distribution of 
given by the concentric fin modeL The moving fin model overpredicts the 
radial heat transfer to the charge in the hot zone while underpredicting 
this heat transfer in the gradient zone. The overpredicted heat transfer 
in the hot zone contributes to <p^ while the lack of heat transfer in 
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the gradient zone reduces ; if the effects of the two errors in 

approximately balance, the moving fin model provides a good 
estimation of . The balance of the two errors is illustrated by the 

dotted curve in Fig. 7.8 which is plotted by reflecting the secticxi of the 
curve for labeled "a" about the axial location C, ■ -.5 ■ - X^/i 

and adding to the section of the same curve labeled "b". If the 
distribution of cedculated by the moving fin model was to exactly 

match the dotted curve, the value of from eq. [7.14a] would 

equal that given by the concentric fin model. Figure 7.8 shows that these 
curves, while not exactly coinciding, are nevertheless close. Therefore, 
the use of the moving fin model approximation for be expected 

to provide a reasonable estimate of at the end of an adiabatic gradient 
zone. Table 7.1 indicates that this conclusion is valid at least over the 
range of K and S examined. 


Xa RftDIAL 






MATURE GRADIENTS GENERATED A T THE GRCWTH mTERFACE IN THE 

Q E .A OSICISLE 


The axial temperature gradient in the charge, , must undergo an 
abrupt change at the growth interface when Pes% ^ 0 or when % ^ 1. (c.f., 
eq. [6.7a]). The axial gradient in the crucible, on the other hand, does 
not change at the growth interface. Figure 7.9 shows that such conditions 
produce a difference between the charge and crucible temperatures at the 
interface which diminishes with distance from the interface. The exchange 
of heat between the charge and the crucible which results from this 
temperature difference has a profound effect on the distribution of in 
the vicinity of the interface. 

Figure 7.10 shows the variation cf in an adiabatic gradient zone 
for systems which are symmetric except Rr ^ 1. The presence of a crucible 
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Fiq. 7.9: Crucible and mean charge ta^ratvire distributions 

interface when Rr 5^ 1 and/or PesRn 0- 
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Fig. 7^ 0; Radial tanperatare variations generated at the 
interf when F 1. 
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is indicated by a value of 1.5 for the diameter ratio S . As % increases, 
it is seen that the distribution of , and more importantly, the value 
of 4^^ at the interface, becomes more negative, causing isotherm shapes 
to be curved in the adverse direction (concave toward the melt). Even a 
small difference in melt and crystal thermal conductivities (e.g., Rr ■ 1.1 
in Fig. 7.10) has a pronounced effect on • 

Values of Rr greater than unity, being typical of semiconductors, have 
been considered in Fig. 7.10. The variation of <p^ which occurs when % is 
less than unity is similar to that shown in Fig. 7.10 except that the 

i 

curves would be displaced toward more positive values of compared to 
the symmetric case of Rr « 1. 

The effect of the generation of latent heat at the growth interface cai 
in the gradient zone is shown in Fig. 7.11. The systems considered are 
symmetric except that PesR^ ^ 0. It is seen that as Rh increases, 
becomes more negative, similar to the effect of Rr > 1 shown in Fig. 7.10. 

Although the curves of Fig. 7.11 are plotted with Rh as a parameter, 
it must be noted that the effect of the generation of latent heat is 
proportional to the product PesRH* The curves of Fig. 7.11 may also be 
interpreted as indicating the effects of changing Pe with a constant value 
of Rfl. The value of Pes » 0.05 used in Fig. 7.11 has been chosen so that 
PesRH would be sufficiently large to have an effect on for typical 
values of Rh. 

The generation of radial gradients at the interface is due to the 
presence of a crucible. Figure 7.12 shows the effect of the crucible 
diameter ratio, S / on for systems which are symmetric except 

Rr = 2. Accordingly, very thin crucibles are required to eliminate the 


hot end 


C in gradient zone 


cold end 


Fig. 7,11 ; The effect of the generation of latent heat 
the radial temperature variations in the oradient zone. 
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Fig. 7.12 ; The effect of the thickness of the crucible on 
the radial tanperature variations generated at the 
interface . 
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interface effect. 

Figure 7.13 shows the effect of changing the location of the interface 
on the distribution of in the adiabatic gradient zone for systems which 
are symmetric except that ■ 2. As the interface approaches the hot 
zone, the value of <p^ at the interface becomes more positive due to the 
influence of radial temperature variations created i.i the hot zone (see 
section 7.2). For the systems considered in Fig. 7.13, however, the 
dominant factor is clearly the conductivity change at the interface. In 
general, the importance of radial temperature variations generated at the 
interface depends on their relative magnitude compared to radial 
temperature variations created by other effects. 

I>iJ. The Two-Fin Thermal Model 

In order to obtain from the results of section 7.1 a useful analytical 
relation which expresses the primary effects of the thermal parameters on 
the radial temperature variatiais near the growth interface, a simple and 
accurate analytical expression for is necessary. Radial heat 

exchange between the charge and crucible in the vicinity of the interface 
is not accounted for in the one-dimensional moving fin model; approximating 
by t-herefore, is not accurate near the interface when 

a crucible is present. The concentric fin model is too complex to yield a 
simple analytical solution for . The present section describes a 

thermal model of the heat transfer '•ear the interface which provides a 
simple approximation for 

The new thermal model is equivalemt to the concentric fin model 
described in (Chapter 6 except for the following: 


(1) The charge is modeled as a single fin as shown in Fig. 7.14. 


hot end 


L, in gradient zone 


cold end 


Fig. 7.13 ; The effect of the interface position on radial 
temoerature v*ariations Generated at the interface. 
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Since there are only two concentric fins — the charge and 
the crucible — the present thermal model is called tiie "two- 
fin" thermal model. The representative radial location for 
the charge temperature is taken as Din/2, defined in Fig. 

6 . 2 , since this location divides the charge equally with 
regard to its cross sectional area. The representative 

radial location for the crucible temperature is again D 3/2 
(c.f.. Fig 6.2). 

(2) Hit' Peclet number is assumed negligible. 

(3) The gradient zone is assumed adiabatic; h«once, tne gradient 
zone annulus is neglected. 

(4) Only the region in the vicinity of the interface ib modeled. 

This region ext«ids from the interface, in either direction, 
at least the distance required for the interface effect to 
significantly die out. The effect of the remainder of the 
growth system is included only through imposing the average 
charge/crucible axial temperature gradients in the gradient 
zone, Gl and Gs in the solid and liquid respectively. These 
parameters are obtained from the one-dimensional moving fin 
model. Radial temperature gradients generated by other 
factors are neglected. 

In accordance with the assumptions listed above, the fin equations ror 
the charge and crucible are: 





■¥ 


t »»\ ~ Tc ir 

^»v»Crj\oc. 


= 0 


[7.17] 
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[7.18] 


The thermal resistance, Rmcr.loc «hioh defines the thermal coupling 
between the crucible and the charge# is given by [30] : 




i\i/x 


[7.19] 


Define the following non-dimensioTal tenperatures; 




'cr 


'cr 




[7.20a] 


[7.20b] 


$ is the average charge/crucible temperature weighted by their respective 
axial thermal condixtances. Using eqs. [7.19] and [7.20] in eqs. [7.17] 
and [7.18] , and nondimensionalizing# results in the followiiig form for the 

twO"fin model fin equations: 


i!£. 


- 0 


[7.21] 



[7.22] 


[7.231 
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Equation [7.21] indicates that <p varies linearly in the gradient zone 
under the present assumptions. This conclusion is also obtained in section 
4.4.3 where the axial gradient of ^ is adso denoted by Gl and Gg in the 
liquid and solid respectively. Section 4.4.3 edso explains that the axial 
temperature gradients in the gradient zone predicted by the moving fin 
model should be interpreted as Gl and Gg. The present two-fin model 
therefore assumes that Gl and Gg are known a priori from the results of the 
moving fin modeL 

The solution of eq. [7.22] is: 


in the liquid: C ^ U] 
in the solid: ^ C,' 

iC-COl [7.24b] 

where the coefficients denoted by C are the unknown constants of 
integration. 

The assumptions of the two-fin model require that 0 as 

; therefore, Cil = ^25 = 0* Continuity of temperature at the 
interface requires ; therefore C 2 L = ^IS- Finally, 

continuity of flux in the crucible fin at the interface requires 

; this boundary condition yields: 

r - f - [7.25] 

^IS " Ws ^ KkVu 
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The difference between Gl and Gg is found from eq. [4.25] : 


UKsCS’--!) 


[7.26] 


Substituting these results into eqs. [7.24] produces the solution for 


cr\ = C d)(pl- Ws (.5- 1; )\ 


[7.27a] 

[7.27b] 


where: 


W& \-^K4[S^-0 


[7.28] 


The constant C defined by eo. [7.28] includes in the numerator the 
superposed effects of the generation of latent heat and the change in 
charge thermal conductivity at the interface. The axial gradient in the 
liquid, Gl, is always negative; if % > 1, as is typical of semiconductors, 
these two effects reinforce each other, causing a larger disturbance at the 
interface than would occur by either effect acting alone. On the other 
hand, if Rr < 1, it may be possible to eliminate the interface effect by 
choosing a value for the product PegRH that would make the numerator of eq. 
[7.28] equal to zero. This could be accomplished by changing either the 
lowering rate (thus changing Peg) or the temperature difference between the 
hot and cold zone furnaces (thus changing Rg). Either possibility, 
however, has potentially harmful side effects on the constitutional 
supercooling requirement. 
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7.3»2 Radial Temperature Variations Near the Growth Interface 

Substituting eqs. [7.27] into the nondimensional form of eq. [7.18] 
yields an expression for due to the heat exchange between the 
charge and the crucible near the interface: 








1 Ki(s‘-0 


C e>cp[-'Ws(,£;-5()] 


[7.29a] 


[7.29b] 


Using eqs, [7.29] in eqs. [7.11] and [7.12] produces the value of at the 
interface due to the interface effect: 


gi, \ C Ws / t»3) ^ Wu / <»2 >L 


[7.30] 


The value of 4*^ at the interface determined by eq, [7.30] is compared 
to the results shown previously in Figs. 7.10 to 7.13 in Table 7.2. In 
general, the comparison is favorable, indicating that the two-fin model 
provides a reasonable approximaticn for in the vicinity of the 
ijiterface. The largest differences occur for those cases which consider 
effects not included in the two-fin model. For example, the cases 
considered in Fig. 7.11 use a larger Peclet number, 0.05, compared to that 
used for the other figures, 0.001, in order to accentuate the liberation of 
latent heat. The effect of the larger Peclet number, not included in the 
two-fin moG?l, is to create a more negative value of • 

A difference between eq. [7.30] and the results from the concentric fin 
model is also apparent for the cases of small S considered in Fig. 7.12. 
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— 1 

eq. [7.30] plus 


variable 

concentric 


attenuation effect 

Figure 

parameter 

. 

fin model 

eq. [7.301 

from eq. [7.9] , 
[7.10], and [7.14] 

7.10 

% - 1.0 

.000000 

.000000 

.000000 


1.1 

-.000493 

-.000464 

-.000464 


1.25 

-.00110 

-.00114 

-.00114 


1.5 

-.00186 

-.00176 

-.00175 


2.0 

-.00274 

-.00261 

-.00260 


4.0 

-.00343 

-.00330 

-.00329 

7.11 

Ru " 0.0 

-.000023 


.000000 


0.25 

-.000219 

-.000179 

-.000179 


0.5 

-.000416 

-.000358 

-.000360 


1.0 

-.000810 

-.000716 

-.000720. 


2.0 

-.00160 

-.00143 

-.00144 

7.12 


-.000021 

-.000002 

-.000021 


1.001 

-.000040 

-.000021 

-.000039 


1.01 

-.000214 

-.000183 

-.000201 


1.05 

-.000818 

-.000730 

-.000744 


1.1 

-.00136 

-.00122 

-.00123 


1.25 

-.00227 

-.00207 

-.00207 



-.00274 

-.00261 

-.00260 

7.13 

H 

o 

o 

-.00274 




-0.1 

-.00245 




-0.2 

-.00210 




-0.3 

-.00166 

-.00252 



-0.4 



-.00108 

-.00249 



The values of Gl and 9\, required for use in 

eqs, [7.14] and [7.30] are listed in Tables 7.3 and 7.4. 


Table 7.2; Comparison of the approximate methods of calculating 

to results of the concentric fin model for the cases 
considered in Figs. 7.10 to 7.13. 
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This error is attributable to the neglect by the two-fin model of radial 
gradients generated in the furnace zones. In order to include an 
approximation of the attenuation effect, the value of (<^i) determined by 
eqs. [7.9], [7.101, and [7.14] is simply added to that determined by eq. 
[7.301. The results of this approximation, shown in the last column of 
Table 7.2, are seen to improve the estimation of for small values of 
S . As S increases, the contribution to due to the attenuation 
effect becomes less important than the radial gradients generated at the 
interface. 

Table 7.2 also demonstrates that eq. [7.301 is unable to account for 
the effects of the changing interface location on This error is 
again attributable to the relative importance at the interface of radial 
gradients generated in the hot and cold zones compared to those generated 
at the interface. As the interface approaches the end of the gradient 
zone, the influence of radial gradients generated in the adjacent furnace 
zone becomes increasingly important. 

The last column of Table 7.2 also provides an approximation for the 
effect of changing interface location on . This correction to 
eq. 17.301 has been calculated assuming that the radial gradients generated 
in the hot and cold zones attenuate at a rate proportional to exp(-7.63C ). 
The actual rate of attenuation is, however, smaller when a crucible is 
present (c.f., section 7.2). It is expected that use of the actual rate of 
attenuation would produce a better correction to eq. [7.301 for the effect 
of changing interface location. Unfortunately, a simple analytical 
expression for the attenuation of from the ends of the gradient zone 
in the presence of a crucible is not available. 
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Cfl c rect i on to the Axial Temperature Gradien t In the Liquid At The 
Growth Tnterfari* 

The moving fin model predicts the value of Gl for the 

axial temperature gradient in the liquid at the interface. The 
constitutional supercooling requirements however, depends on the value of 
at the interface. The difference between these temperature 
gradients can be approximated from the results of the two-fin model. 
Differaitiating eqs. 17.20] with respect to ^ and solving for 
yields: 






1l '■ 

d? ] 


[7.31] 


Differentiating eq. [7.27a] with respect to Q and substituting into 
eq. [7.31] yields: 


'UAZi ) 

AC 


- Qu 




w^C 


[7.32] 


The second term on the right side of eq. [7.32] is a correction to be 
applied to Gl in order to account for the effect of radial heat transfer 
near the interface on the axial gradient in the charge at the interface. 

The axial temperature gradient in the liquid at the interface 
calculated by the moving fin model, Gl, and by eq. [7.32] are compared to 
results from the eccentric fin model in Table 7.3 for the cases considered 
in Figs. 7.10 to 7.13. It is seen that eq. [7.32] provides an excellent 
approximation to even though the correction term may be 
significant (e.g., large % or large R{j). The approximation is most in 
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axial temperature gradient in the 
liquid at the growth interface 


variable 

concentric 

Figure parameter 

fin model 


model 

Gl 


-.350 

-.340 

-.326 

-.306 

-.272 

-.190 


-.350 

-.332 

-.307 

-.275 

-.228 

-.140 



Table 7.3: Comparison of the approximate methods of calculating 

to results of the concentric fin model for the cases 
considered in Figs. 7.10 to 7.13. 
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error as the interface approaches the hot zone (c.f., C; ■ 

Fig. 7.13) due to the additional influence of radial heat trar.crer in the 

hot zone. 

7.3.4 fprror^inn Tn th«> Interface Temoeratura 

The moving fin model predicts for the interface temperature the 
average charge/crucible temperature* i.e.* &L interpreted as • 

•Rie difference betwe«i ^ (Ci) 3nd ^ (ijj) found by the two-fin model may be 
used as a simple correction to the value of found by the moving fin 

model: 

Bi - 4>r. ( 5.0 - (50 - 

Rearranging eqs. [7.20] yields: 







- Ks (&^--0 

\+KsCS’--0 



[7.34a] 


[7.,'>*.b] 


Since both and have the same value in the liquid and solid 

at the interface* eqs. [7.34] imply that ^ 

two-fin model predicts a step ch 2 uige in at the interface. This is due 
to the definition of ^ which averages the charge and crucible 
temperatures according to their respective axial thermal conductances. 
Since the moving fin model neglects the temperature difference between <p^ 
and it does not predict this step change in at the interface. 

In order to approximate the effect of radial heat transfer near the 
interface on the interface temperature* the correction term within brackets 
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in eq. C7.33] is taken simply as the arithmetic average of eqs. [7.34a] and 
[7.34b]. Using eq. [7.27] evaluated at 5 yields: 


- 0 ;, + 





K,(S^-0 1 

'^i^fKsCS'-O i 


[7.35] 


Uie interface temperature determined by eq. [7.35] is compared to that 
found from the concentric and moving fin models in Table 7.4 for the cases 
considered in Figs. 7.10 to 7.13. It is seen that the moving fin models 
provides a reasonable estimation to the interface temperature, , 
without the correction afforded by eq. [7.35]; this indicates that radial 
heat transfer at the interface between the charge and the crucible has a 
relatively small effect on the interface temperature or location. 
Utilizing the corresponding values of Gl listed in Table 7.3, the moving 
fin model would mislocate the interface by at most 0.1 charge diameter for 
the cases considered (c.f., Rr » 4 from Fig. 7.10). Results calculated 
from eq. [7.35], however, compare yet more favorably to those of the 
concentric fin model. 

The largest errors in the estimate of the interface temperature 
considered by Table 7.4 are attributable to the larger Peclet number used 
in the cases of Fig. 7.11. Calculating through the use of the non~ 
zero Peclet number relation, eq. [B.8] , increases the accuracy of 
eq. [7.35]; for the cases considered by Fig 7.11, the difference 
between 9 ^ calculated in this way and by the concentric fin model is no 
more than one digit in- the third significant figure. 
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teirperature of *±e growth interface 


variable 

concentric 

moving fin 


Figure 

parameter 

fin model 

model 

eq. [7.35] 




(Pe » 0) 


7.10 

Rr • 1.0 

.500 

.500 

.500 


1.1 

.509 

.507 

.509 

• 

1.25 

.522 

.518 

.522 


1.5 

.540 

.534 

.540 


2.0 

.570 

.561 

.572 


4.0 

.648 

.634 

.651 

7.11 

Rg “ 0.0 

.517 

.500 

.500 


0.25 

.522 

.504 

.505 


0.5 

.527 

.508 

.509 


1.0 

.536 

.516 

.519 


2.0 

.554 

.532 

.537 

7.12 

% - 1.0001 

.619 

.619 

.619 


1.001 

.619 

.618 

.618 


1.01 

.617 

.617 

.617 


1.05 

.610 

.610 

.611 


1.1 

.603 

.602 

.604 


1.25 

.587 

.584 

.588 


1.5 

.570 

.561 

.572 

7.13 

» 0.0 

,570 

.561 

.572 


- 0.1 

.602 

.593 

.604 


- 0,2 

.633 

.624 

.635 


-0,3 

.664 

.655 

.665 


-0.4 

.693 

.685 

.695 




The values of Gl required for use in eqs. [7.35] are 
listed in Table 7.3. 


Table 7.4; Comparison of the approximate methods of calculating 6^ 


to results of the concentric fin model for the cases 
considered in Figs. 7.10 to 7.13. 
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2^ DTABATIC GRADIF VT ZQlffi ^ 

Two additional parameters are required by the concentric fin model 
when heat transfer in the gradient zone is considered: Bi^ and R<;. The 

Biot number in the gradient zone, Bi(;, indicates the thermal coupling 
between the gradient zone annulus and the charge. R<;; is the ratio of the 
axial heat conducting capacities of the gradient zone annulus and the 


charge (c.f., eqs. [6.21). The one-dimensional model of a diabatic 
gradient zone developed in section 4.4.4 shows that the gradient zone 
approaches abiadatic conditions as either Bi(; or Rq approaches zero. 
Otherwise, the temperature difference between the gradient zone annulus and 
the charge ( ^ produces heat transfer which affects axial 
temperature gradients and also, as shown in this section, radial 
temperature gradients within the charge. 


7.4.1 Results From the Concentric Fin Model 

Figures 7.15a, 7.15b, and 7.15c show the variation of in a 
diabatic gradient zone as calculated by the concentric fin model. The 
systems in these figures are symmetric and do not consider a crucible. 
Since the systems are symmetric, the curves have been plotted only for the 
hot side of the gradient zone. Each figure considers a constant Biot 
number throughout the furnace while the separate curves in each figure 
consider various values of R< 5 . The graphs inserted in each figure are 
axial temperature distributions of and corresponding to the 
individual curv® of . As a reference, the distribution of \ in an 
adiabatic gradient zone (i.e., Bi^ = 0) is also included in each figure. 

According to Figs. 7.15, the value of P <5 significantly affects the 
distribution of in the gradient zone. When Rq is small, heat exchange 
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between the charge and the gradient zone annulus occurs only near the ends 
of the gradient zone; toward the center of the gradient zone is nearly 
equal to . In such cases, the distribution of <p^ is only slightly 
altered from that in the corresponding system with a.* adiabatic gradient 
zone. As Rq increases, tlie temperature difference 4?^^ - <p^ is more 
pronounced and increases the radial heat transfer to the charge; -"^s a 
consequence, <p^ increases as welL For sufficiently large values of Rq, 
is approximately linear; further increases in Rq do not affect 
either - c£i or d> but serve only to conduct more heat from the hot to 
the cold zone furnaces through the gradient zone annulus. 

Figure 7.16 demonstrates the effect of increasing BIq for systems 
which have a linear variation of <p^ c- calculated by the concentric fin 
modeL A linear variation of <p^ has been obtained by choosing a large 
value of Rq (i.e., R(2=100i,. Increasing Bi^ has a slight effect on the 
axial temperature distribution in the gradient zone, tending to reduce 
driving force for radial heat transfer in the gradient zone, <p^ - . 

This effect, however, is more than offset by the augmented thermal coupling 
between the charge and gradient zone annulus represented by the larger 
value of BiQ. As shown in Fig. 7.16, the result is that increases with 
BiQ. The curve for Bic = 5 shows, for example, that <p^ increases, rather 
than attenuates, from the end of the gradient zone when the gradient zone 
Biot number is greater than that in the furnace zones. 

Figure 7.17 shows that the effect of heat transfer in the gradient 
zone when a crucible is present is similar to those systems examined above 
without a crucible. (The system conside.red in Fig. 7.17 is analogous to 
that of Fig. 7.15b except for the presence of a crucible.) 
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7.4.2 anaiviiical /xnprnvimfcinn for the Qt Hprit Tr a nafftc In . t to. 

Zone 

The effect of radial heat exchange between the furnace and the cha.rge 
on the value of at a parUculac axial location is described by 
eq. (7.41. Analytical expressions for the effects of radial tebperature 
gradients created in the hot and cold zones and at the growth interface 
have been developed using eg. 17.41 in sections 7.2 and 7.3 respectively. 
The present section considers a simple model of heat transfer to or from 
the charge in the gradient zone in order to analytically approximate its 

effect on . 

The present model considers thr.“ the primary effect of heat transfer 
between the gradient zone annului t e charge is attributable to thermal 
conditions neat the axial looaao-. -r interest. This assumpUon is 
qualitatively justified by the high rate of attenuation of the effects of 
radial heat transfer distant from the locaUon of interest: c.f., 

eq. (7.41. Accordingly, a lineatized description for 
location of interest, C f is employed: 

^ r) tV.36] 

^ <ix 1. 

Substituting eq. [7.361 into eq. [7.71 yields: 

, \ [7.371 

Next, eq. [7.161 is used to approximate > this approximation 

assumes that of the moving fin model accurately rf^presents 

of the concentric fLj model and that i-Tjc 
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linearly near the locatioi of interest. Subsituting eg. [7.16] into [7.37] 
yields: 

(P (z) = [7.38] 

14 . 55 ^ 

(WithLn the gradient zone, <P^ becomes .) According to eg. [7.38] , the 
value of at a certain location, C , depends only on the local Biot 
number and the local temperature difference between the furnace and the 
charge. 

Eguation [7.38] is expected to be a good approximation of <t>^ when 
is approximately linear in the vicinity of the axial location of 
interest. Linearity of is not expected, however, near the ends of 

the gradient zone (see section 7.2) or at the interface (see section 7.3). 
In these regions, eg. [7.381 should be interpreted as an order of magnitude 
estimation of the effect of heat exchange in the gradient zone on the value 
of which would exist in the absence of such heat exchange (i.e., due to 
the attenuation and interface effects) . 

7.4.3 Constant Inter !race Sh a pe in the Grad i ent zone 

Eguation [7.38] suggests, in accordance with its limitations discussed 
above, that may be maintained at a ccaistant value over a portion of the 
gradient zone if the product BiQ*( ~ ^ constant over a 

sufficiently long region near the center of the gradient zone, i.e., not 
near the ends of the gradient zone. The following set of "design rules" 
illustrates how the heat transfer in the gradient zone can be "tailored" so 
that these ccnditions may be approximately achieved for a system in which 
the interface di.sturbance is not present: 
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(1) Calculate in the gradient zone, assuming that Bi <3 ■ 0, by 

using the results of the moving fin model. A symmetric system is 
presently used to demonstrate this procedure; consequenty, eq. [5.8] 
provides the required value of • 

(2) Allow Rq to be sufficiently large so that (p^ is linear. 

(3) Alter the end temperatures of the gradient zone annulus so that 

. In this way, 4l~<P is constant in the gradient 
zoiie, negecting the effects of the crucible near the aids of the gradient 
zone. In the present development, it is arbitrarily chosen to maintain 
while altering • In thij case: 






[7.39] 

[7.40] 


(4) The desired value of determines the required value of Bi^ by 
substituting eq. [7.40] into eq. [7.38]: 




[7.41] 


A non-zero value of BIq* leads to a non-symmetric system and non-constant 
( )f thu' contrc’dicting the assumptions of steps (1) and (3) above. 
These effects are shown, however, to be of secondary importance in the 
results which follow. 

(5) The desired radius of curvature of the interface, N^, is next 
used to eliminate <p^ from eq. [7.41]. Using eq. [G.6] ■ 



2^.1 


original 

OF POOR 


Page 

Quality 


[7.42] 


Using eq. [5.8] to eliminate 


;»• = - i2lL 17.431 

Hiree teSw cases are employed to demonstrate the procedure described 
above. The parameters defining the systems and the results of the 
calculations of steps (1) through (5) are -shown in Table 7.5. 
Figures 7.18, 7.19, and 7.20 plot the resultant distributions of 
temperature and isotherm curvature calculated by the concentric fin model: 
part (a) of each figure plots the distribution of in the gradient zone 
while part (b) shows the distribution of <fe and S and the inverse 
radius of curvature, N"l, of isotherms in the gradient zcxie. The desired 
constant radius of curvature chosen for use in eq. [7.43] is Ni = 20; 
therefore 0.05. 

Figure 7.18 considers a system without a crucible and with a gradient 
zone length = 1.0. It is seen that the distribution of and N tend 
to flatten out near the center of the gradient zone, approximating the 
desired values of » 0.00127 and 0.05. It is expected that the 
distributions of and N will be yet more flat when the influance from 
the ends of the gradient zone are more diminished at its center. This may 
he accomplished by increasing the length of the gradient zone to A^ “ 2.0 
(Fig. 7.19). The curves for and N are quite flat in a large portion of 
the center of the gradient zone in this case. Note that the required value 
of <p^ , 0.000903, for case 2 (see Table 7.5) is larger than actually 
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Case 

1 

2 

3 

Figure 

7.18 

7.19 

7.20 

S 

1.0 

1.0 

1.5 


1.0 

1.0 

2.0 

Big » Bi<; 

0.5 

0.5 


Big* » Bic* 

0.471 

0.471 

0.291 

ae«/jc 

(eq. [5.81) 

-0.407 

-0.289 

-0.259 

(eq. [7.39]) 

0.593 

0.422 

0.481 

- 4>m) 

(eq. [7.40]) 

0.297 

0.211 

0.241 

_ . ** 

1 Big 

(eq. (7. 43]) 

0.0625 

0.0625 

0.0491 

BiG 

0.0625 

0.0625 

0.0735 


0.00127 

0.000903 

0.000811 


Synmetric system except ^ 0 

“ Cod 

Desired value of is 20 

Table 7.5; Description of the systems used in Figures 7.18r 
7.19, and 7.20. 
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Bi * 0.5 
BIq* 0.0625 
No crucible 

Symmetric system except- 
= 0.593 


hot end 


C in the gradient 
zone 


cold end 


'^•I8a ; Demonstration of the procedure used to obtain 
an approximately constant isotherm shape in the gradient 
zone. Jfo crucible, syirmetric system and \c, = 1.0. 

Distribution of p,-. . 





3iu * B'r * 

page rs 

n w 

BIq = 0.0625 

Oi POOR QUALITY 

No crucible 


Symmetric system except: 


<|)q(Xq/2) * 0.593 



r in the gradient 
zone 


cold end 


Fiq. 7.18b: DeitiDns'cration of the procedure used to obtain 

an approximately constant isotherm shape ^ the gradient 
zone'. NO crucible, symmetric system and = 1.0. 

Distribution of 4>m/ ^tid N. 






F ig. '^.19a ; Analogous to Figure 7.18a except the 'pr'^dient 
zone is longer; i.e., \r; = 2.0. 

Distribution of o • . 


- 192 - 




7one 

us to Figure 7.18b except the 
ger ; i.e. , = 2 . 0 . 

5(3 and N. 







Fig. 7.20a : Analogou?,- to Fic.rure 7.19a except a crucible 

is present; i.e., 5 " 1-5. 

Distribution of <t>A. 
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present as shown in Fig. 7.19. This is due to heat transfer to the charge 
in the gradient zone, not considered in eq. [7.40], which causes <p^ - 
to be slightly smaller than the value predicted for a symmetric system. 
Reduced heat transfer to the charge, and therefore a reduced value of 
results. The value of can be adjusted upward by an appropriate 
alteration in Bi<;. 

Figure 7.20, which considers a system similar to that of Fig. 7.19 
except S * 1*5, demonstrates that the presence of a crucible is adequately 
accounted for in this procedure. 

The procedure outlined i this section is not intended as a 
generalized design method leading to optimized heat transfer for Bridgman 
crystal growth. Systems of pra ical interest do not behave as those 
considered in Figs. 7.18 to 7.20 due to the generation of radial 
temperature variations at the growth interface (c.f., section 7.3). Its 
intent, rather is to accentuate three important features concerning the 
radial temperature distribution in the charge; 


(1) Heat transfer in the gradient zone is a useful tool for 
altering the interface shape in the gradient zone. 

(2) ^There appropriate, results from one-dimensional models, in 
conjunction with the Green's functiai approach, can be used 
to help predict radial temperature variations in the charge. 
In the model developed in this section, eq. [7.38] is seen to 
provide at least order of magnitude accuracy. 

(3) Generation of radial temperature variations at the interface 
(Rr # 1 and/or PesRn 0^ is detrimental for achieving a 
region about the interface which exhibits a relatively 
constant isotherm curvature. 
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a.l AXIAL TEMPERAmjRE DISTRIBUTION 

The principal objecvives of the thennal design of a Bridgman growth 
system regarding the axial temperature distribution of the charge are 
(1) to obtain a sufficiently large axial temperature gradient in the melt 
at the growth interface and (2) to control the location of the interface. 
The thermal parameters affecting the attainment of these objectives may be 
classified as fixed/ i.e./ the thermal properties of the charger or 
variable/ i.e./ the furnace process parameters such as the lowering rate/ 
diameter of the charge/ diameter and thermal conductivity of the crucible/ 
gradient zone length/ furnace heat transfer coefficients and furnace 
teiperatures. 


a>l.I Thermal gcopertiea q£ the Qarqe 

The thermal properties of the charge are fixed parameters and 
therefore constrain the thermal behavior of the growth system. The most 
important of the thermal properties of the charge is its thermal 
ccxiductivity. A small thermal conductivity leads to a large Biot number; a 
Bridgman system for the growth of CdTe, for example, can have a Biot number 
in excess of unity. In such caseS/ satisfying the constitutional 
supercooling requirement may not be difficult. Additionally, since the 
"infinite" charge length is smaller when the Biot number is larger, a 
charge with a small thermal conductivity exhibits decreased thermal 
transients associated with growth near the ends of the charge and therefore 
tends to stabilize the growth interface. Biot numbers are smaller, on the 
other hand, when the thermal conductivity of the charge is relatively large 
(e.g., germanium). In this case, it is important to verify that the axial 
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temperature gradient at the growth interface is sufficiently large; if not, 
efforts to increase the Biot number by augmenting the heat transfer 
coefficients may be necessary. Further, a longer charge and furnace is 
necessary for the growth of an equivalent length of crystal under quasi- 
steady thermal conditions when the thermal conductivity is relatively 
large. 

Related to the magnitude of the thermal conductivity is the change in 
thermal conductivity between the melt and the crystal at the growth 
interface (denoted by the symbol Rr). When the melt conductivity is 
greater than that of the crystal (i.e., Rr > 1) the axial gradient in the 
liquid at the interface is adversely affected, becoming less than that on 
the solid side of the interface. The magnitude of this effect may be 
significant (c.f., section 5.3) and is a potential problem for the growth 
of semiconductors which typically have values of Rr greater than unity. 
When Rr is less than unity, as is the case for metals, satisfying the 
constitutional supercooling requirement is more easily effected as a 
consequence of this effect. 

Other thermal properties of the charge which affect its axial 
temperature distribution are the thermal diffusivity, , mass density, 
^ , and latent heat of solidification, AHg]_. The thermal diffusivity is 
a factor in the Peclet number which accounts for the effect of convection 
of sensible heat due to the motion of the charge. The product of ^ and 
AHgi is contained in the product of nondimensional parameters PesRfl which 
accounts for the effect of the liberation of latent heat at the growth 
interface. Additionally, the convection and latent heat effects are 
proportional to the lowering rate, V. The significance of these heat 
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flows, however, depends upon their relative magnitude compared to the 
primary heat flow in the system, i.e., the heat flow between the furnace 
and the charge in the hot and cold zones. These comparisons are made in 
eqs. [5.6] and [5.121 where indicates the importance of the furnace- 
to-charge heat transfer. The effects of the thermal properties ^ , f , and 
AHsi therefore depend on the prevailing values of Bi and V; typical values 
of these parameters for semiconductor growth systems produce small 
convection and latent heat effects. It may be concluded, therefore, that 
these thermal properties are of secondary importance to the axial 
temperature distribution of the charge compared to the thermal condix:tivity 
of the charge. 

AJLa 2 Furnace giogesa. Pacameteca 

A number of furnace design parameters can be used to alter the axial 
temperature distributicxi of the charge. The most obvious of these are the 
hot and cold zone furnace temperatures. An increase in the hot and cold 
zone furnace temperature difference results in a roughly proportional 
increase of the axial temerature gradient within the charge. On the other 
hand, lowering one or both of the furnace temperatures moves the growth 
interface toward the hot zone while an increase moves the growth interface 
toward the cold zone. Furthermore, alterations of the axial gradient and 
interface location can be made independently of each other by the proper 
choice of furnace temperatures; for example, the axial gradients can be 
made to increase while maintaining the interface in the same location or 
vice-versa. Typically, temperature limits in the hot and cold zone furnace 
constrain the application of these techniques. The heat pipes used in the 
experimental system described in section 2.2.1 may bo used only between the 
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temperatures of 500 *C and 1100 *C. Furthermore, the hot zone furnace 
temperature may be limited by the development of excessive vapor pressure 
of one of the constituent elements in the melt or by the mechanical 
properties of the crucible at elevated temperatures. 

The most important of the dimeisiailess parameters regarding the axial 
temperature distribution of the charge are the Biot numbers of the hot and 
cold zones. The Biot numbers reflect the thermal coupling between the 
furnace and the charge; eq. [5.91 indicates that as the Biot numbers 
increase, so do the axial temperature gradients in the gradient zone. The 
effect of the thermal caiductivity of the charge on Bi has beai discussed 
in section 8.1.1. The furnace process parameters which affect the Biot 
numbers are the furnace to charge heat transfer coefficients, h, charge 
diameter, D, and the thermal conductivity emd thickness of the crucible, K 
and S through their effect on the effective Biot number, Bi*. It must be 
noted, however, that the internal radial thermal resistance of the charge, 
accounted for by Bi**, limits the degree to which axial gradients nay be 
increased through augmeitati'n of the Biot number. As shown by eq. [4.14], 
Bi** can not exceed the valued of 8; changes in h, K and & have little 
effect on increasing Bi** when it is close to this maximum value. 

The value of the furnace to charge heat transfer coefficient, h, 
includes both radiative and conductive components. The radiative component 
may be minimized by employing reflecting surfaces between the inner furnace 
and outer crucible surfaces; it may be increased by utilizing surfaces of 
large emissivity. The magnitude of radiation heat transfer is, however, 
limited to that which would occur between black surfaces. Additional 
amounts of heat transfer may be obtained by conduction heat transfer across 
the furnace cavity gap. Appendix C indicates that the conductive part of 


- 200 - 


QiH 


rP'C'v-. r-. . 

cr ■ OrV Qi/AUTY 

the gap heat transfer can be as large as the radiative part by using an 
interjacent gas of high thermal conductivity^ such as helium, and a gap 
width of the order of one to three millimeters. Higher conduction heat 
transfer coefficients may, of course, be obtained by further decreasing the 
gap width. Small gap widths accentuate, however, the error in misaligr’^'^nt 
of the axes of the charge and furnace cavity which lead to non-radially 
symmetric heat trSi-iSfer coefficients and interface sh\pe. (The radiation 
component to the gap heat transfer is less sensitive to the gap width, 
especially for gap widths which are small compared to the furnace cavity 
diameter. 

The thermal coupling between the furnace and the charge, is in most 
cases adversely affect«<3 by the crucible; i.e., the effective Biot number, 
Bi* is usually less than the Biot number, Bi (see section 4.2.2). In order 
to maintain a large value of Bi*, it is desirable to use a crucible which 
has a small thickness compared to the charge diameter; i.e., S should be 
as close to unity as possible while maintaining structural integrity. The 
use of crucibles which have a very high or very small thermal ccaiductivity 
compared to that of the charge should also be avoided since they tend to 
cause a significant reduction in Bi*. Equation [4.7] may be used to reveal 
which among several potential crucible materials will yield the largest 
effective Biot number. 

An increase in charge diameter, D, tends to increase the dimensionless 
axial gradient in the gradient zone through its effect on the Biot number. 
As section 5.6 demonstrates, the corresponding effect on the dimensional 
axial gradient depends on the prevailing values of Bi, K and S . Unless 
accompanied by a decrease in S f an increase in D causes a decrease in the 
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axial temperature gradient at the interface. A decrease in <S , however, 
does not guarantee an increased axial gradient; such cases should be 
evaluated individually. 

The length of the gradient zone, , has an appreciable effect on the 
axial temperature gradients at the interface. For example. Fig. 5.3 
demonstrates that the axial gradi^t may always be increased by decreasing 
A^. The effect is most significant when the hot and cold zone Biot 
numbers are large (greater than about 0.1). The absence of a gradient 
zone, (i.e.,. A^^ » 0), however, limits the extent of this effect; if the 
axial temperar.ure gradient is not sufficiently large when A^ « 0, other 
methods to increase the axial gradient must be employed. Hie gradient zone 
length additionally affects the distribution of radial temperature 
varitions near the growth interface. "Changes in A^ must therefore be 
evaluated on the basis of its effect on the interface shape as well as on 
the axial temperature gradient at the interface. 


radj:al 




MATURE VARIATIONS NEAR THE GRCWTO 


RFACE 


Radial temperature variations within the charge have been shown in 
Chapter 7 to attenuate rapidly with axial distance from the location at 
which they are generated. This behavior permits the formulation of a 
thermal design procedure leading to a satisfactory interface shape for 
systems in which radial temperature variations created at the interface are 
not present. In summary; 


The interface should be located in a gradient zone which is 
made approximately adiabatic through the use of a thermal 
insulating material placed between the hot and cold zones and 
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radiation shields in the farnace cavity gap. The gradient 
zone should be sufficiently long so that radial temperature 
gradients generated in the hot and cold zones have largely 
attenuated at the interface location; this length is 
approximately cxie to two charge diameters and depends on the 
effect of the crucible in reducing the rate of attenuation of 
the radial temperature gradients created in the furnace zones 
(see section 7.2.2). The exact location of the interface 
within the gradient zone is chosen so that the shape of the 
interface is slightly curved toward the solid. Additionally, 
small amounts of heat transfer to the charge in the gradient 
zone, as described in sections 7.4.2 and 7.4.3, can be used 
alter the interface shape or to maintain a large axial regi(xi 
of nearly constant isotherm curvature. 

The thermal design is complete at this stage if the axial temperature 
gradient in the melt at the interface is sufficiently large to satisfy the 
constitutional supercooling requirement. If the axial gradient is not 
sufficiently large, however, techniques to increase its value, as described 
in section 8.1, can be evaluated while attempting to maintain satisfactory 
interface shape. 

Radial temperature variations generated at the interface are, however, 
the dominant factor with regard to the interface shape. The procedure 
outlined above, therefore, does not lead to desirable interface shapes 
unless the interface effect can be eliminated. This ; '•>sue is of particular 
importance for the growth of semiconductor materials which have values of 
Rr greater than unity; for such materials, the interface effect tends to 
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produce an interface shape which is curved in the adverse direction, i.e., 
toward the melt. 

Counteracting the Effects of Radial Tartperature Variations Generated 

at the Growth Interface 

The analysis of section 7.1 suggests that heat transfer from the 
furnace to the charge t^ds to produce isotherm shapes which are concave 
toward the solid. Such heat transfer in the gradient zcne thus serves to 
make the interface shape more favorable if it is originally concave toward 
the melt. The gradient zone, in such a case, resembles somewhat an 
extension of the hot zone furnace. The potential for completely 
eliminating an adverse interface curvature may be estimated by comparing 
the value of produced at the interface due to the interface effect 
(eq. [?.30]) to a representative value of caused by heat exchange with 
the furnace at the prevailing Biot numbers. The value of at the end of 
an adiabatic gradient zone (see section 7.2) may be used to estimate this 
latter value of . If radial temperature gradients caused by heat 
transfer in the furnace zones are of comparable magnitude, or larger, than 
those generated due to the interface effect, additional heat transfer to 
the charge in the gradient zone may prove useful in reversing an adverse 
interface curvature. If, on the other hand, radial temperature gradients 
caused by heat transfer to the charge are less than those generated at the 
interface, this technique can not be expected to change the direction of 
the shape of the interface. 

A decrease in the value of £ likewise decreases the magnitude cf the 
radial temperature gradients generated at the interface, as shown in 
Fig, 7.12. The value of <5 can be decreased by either an increase in the 
charge diameter or by a decrease in the thickness of the crucible. A 
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decrease in S , while not eliirlnating the interface effectr may be useful 
in combination with heat transfer to the charge as described in the 
previous paragraph. 

The ideal method of dealing with the interface effect is to eliminate 
it. If this is possible^ the design procedure for obtaining the desired 
interface shapef described in section 8.2, becomes applicable. The cause of 
the interface effect is the change in the axial temperature gradient within 
the charge at the growth interface which is not present in the crucible. 
Efforts to eliminate the Interface effect therefore must either (1) remove 
the change in axial gradient within the charge at the interface or (2) 
produce an equivalent change in axial gradient within the crucible. 

Equation [7.26] suggests that the slope change at the interface 
disappears for em appropriate value of the lowering rate, V, when < 1. 
T!he efficacy of this technique depends, however, on the interplay between 
the value of V required to elir^nate the slope change at the interface and 
tlie effect of V on the constitutional supercooling requirement. Further, 
this technique is not applicable wh^ % > 1. 

Since the electron potentials are not equal in the solid and liquid 
portions of the charge, the passage of an electric current across the 
growth interface produces or absorbs heat (depending on the direction of 
the current) by virtue of the Peltier effect. Such heat can be used to 
eliminate the change in axial temperature gradient at the interface caused 
by the change in thermal conductivity and the liberation of latent heat. 
Rewriting the h'at balance at the interface in order to include the Peltier 
heat yields: 
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where: ■ Gg in order to eliminate the slope change at the 

interface 

Se ■ Seebeck coefficieit 

I* ■ current per unit cross sectional area 

SeT^I” ■ rate of generation of Peltier heat per unit 
cross sectional area 

Tj, ■ absolulte melting temperature 


Using germanium as an example (Se - 70^V/Rr Tj, ■ 1211*K), eq. [8.1] shows 
that the required curreit density is in excess of 200 eunps/cm^ for typical 
operating conditions. Such large current densities may be difficult to 
attain experimentally and may also produce significant Joule heating 
effects. Although each case should be examined individually^ it appears 
that Peltier heating is not a viable method for eliminating the slope 
change at the interface. 

Peltier heating/cooling has, however, greater potential for modifying 
the interface shape in small amounts. For example, a flat interface shape 
may be made slighlty ccxicave toward the solid by the applicatioi of a small 
amount of Peltier cooling. This technique yields a similar result as that 
of transferring small amounts of heat to the charge in the gradient zone as 
described in section 7.4.3. 

The interface effect may also be eliminated by producing a slope 
change in the crucible at the interface which is equivalent to that in the 
charge. This may be done by adding the appropriate amount of heat to the 
crucible in a very narrow region near the axial location of the interface. 
The remainder of the gradient zone may be adiaba.ic if desired. A thin 
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The objectives of this analysis are to determine the primari’ thermal 
behavior of the vertical Bridgman crystal growth method and, where 
possible, to develo^j simplified analytical results to indicate 
quantitatively the relationship between the thermal parameters of the 
system and its thermal response. The purpose of these objectives are (1) 
to aid in both the preliminary thermal design of a Bridgman growth system, 
(2) to better interpret experimoital results of existing systems in order 
to improve their design and (3) in a more general sense, to critically 
evaluate the capabilities of the Bridgman technique. The major efforts 
involved in attaining these objectives has been the development of 
appropriate simplified thermal models of the Bridgman growth system. The 
major findings of this study are reviewed, and possible directions for 
future research are outlined, in this chapter. 

iJ. gpIKLUSlONS 

A one-dimensional thermal model was developed in order to determine 
the axial temperature distribution within the charge. Approximations for 
the radial thermal resistance offered by the crucible and by the charge 
itself are incorporated in a straightforward manner into the one- 
dimensional model through the definition of the "effective" Biot number. 
In general, these radial resistances decrease the effective thermal 
coupling between the furnace and the charge, leading to reduced effective 
Biot numbers and axial temoerature gradients at the growth interface. The 
relationship betwf the crucible thickness and thermal conductivity and 
the effective Biot number is analytically described by eq. [4.7], This 
relati'' \ship demonstrates that the ratio of thermal conductivity of the 
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crucible to that of the charge, K, should be between 0.1 and 1.0 in order 
to prevent severe reduction of the effective Biot number. E)quatic»i [4.14] 
incorporates the effect of the radial thermal resistance within the charge 
and indicates that axial gradients are limited since the effective Biot 
number can not exceed the value of 8, Efforts to increase the heat 
transfer coefficient between the furnace and crucible therefore does not 
increase axial temperature gradients in systems which already have an 
effective Biot number which approaches this value. 

The axial temperature distribution within the charge is easily 
obtained from the oie-dimensicnal model by computer. Analytical solution, 
however, yields several important results. Criteria are developed which 
indicate wh&i the effects of axial convection of sensible heat, generation 
of latent heat and length of the charge can be neglected. Additionally, 
analytical expressions for the axial temperature gradient in the melt at 
the growth int :fa< (required for the constitutional supercooling 
requirement) and the axial location of the interface are developed. (When 
the axial convection of sensible heat is sufficiently small, these 
analytical results reduce into yet more simple forms.) Hiese expressions 
assume a charge length which satisfies the infinite length criteria but 
includes the effects of; hot and cold zone effective Biot numbers, unequal 
melt and crystal thermal conductivity, axial convection of sensible heat, 
generation of latent heat, length of an adiabatic gradient zone, and 
interface position within the gradient zone. The effects of the charge 
diameter and thermal properties are included by virtue of the non- 
dimensional formulation of the thermal model. Such results indicate that 
the principle factors affecting the axial temperature gradient at the 
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interface are; level of the effective Biot numbers, unequal melt and 
crystal thermal conductivity, gradient zone length and generation of latent 
heat. 

Two-dimensional thermal modeling is utilized to predict the shape of 
the growth interface. The two-dimensional thermal models developed here 
are extensions of the one-dimensional model; i.e., the charge is considered 
as two radially concentric fins, the crucible is a third fin and, in the 
gradient zone, the furnace is considered to be a fourth annular fin. This 
approach is taken, rather than a finite difference or finite element 
computer formulation, for the pot«itial to develop analytical results. 

The equations describing the temperature distributions of the 
concentric fins indicates that radial temperature variations within the 
charge are created by heat exchange between the charge and its 
surroundings (see eq. (7.2]}. Heat transfer to the charge leads to 
interface shapes which are concave toward the solid, and heat transfer from 
the charge tends to produce an interface shape which is concave toward the 
liquid. The effect of heat transfer to the charge at a particular axial 
location dies out rapidly with distance from that location. Thermal 
conditions as little as one charge dieuneter from the interface have only a 
very small effect on its shape. The design of a Bridgman apparatus for the 
proper interface shape ntay therefore be interpreted as providing the proper 
amount and distribution of heat exchange to or from the charge in the 
vicinity of the interface. 

The two-dimensional results identify for the first time a potentially 
serious problem for the satisfactory growth of materials, such as 
semiconductors, which have a thermal conductivity of the melt which is 
greater than that of the solid. In the presence of a crucible, inherent to 
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the Bridgman technique, such conditions lead to heat transfer from the 
charge to the crucible. This is the dominant effect on the shape of the 
interface since it occurs at the interface itself. Further, the shape of 
the interface produced through this effect is curved in the undesirable 
direction, i.e, toward the liquid. 




REQDMMENDATIQNS FOR 




RESEARCH 


The results of the one- and two-dimensicaial models presented in this 
work compare well to published analytical results from more accurate 
computer models. Published experimental results are, however, few. and not 
suitable for the direct verification of the conclusions presented here. 
Such verification is considered necessary and is, in part, the motivation 
for the experimental program currently in progress. 

The numerical accuracy of the simple thermal models of this work is 
expected to be sufficient for initial furnace design as well as to indicate 
the expected parametric thermal behavior. Increased numerical accuracy 
obtainable from a two-dimensional finite difference or finite element 
computer model will, however, be required for final system design and for 
quantitative comparison with experimental results. Effects such as thermal 
conductivity change and/or generation of latent heat at the interface in 
the presence of a crucible, presently not included by computer models 
described in the literature, must be included in any new computer model. 
Modeling assumptions used in the present work, but anticipated to have a 
minor effect in actual growth systems, can also be relaxed; these include 
temperature dependent thermal properties and heat transfer coefficients and 
temperature uniformity of the hot and cold zone furnaces. Finally, the 
furnace boundary conditions should be flexible, especially in the gradient 
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zone, in order to study the effect of different distributions of heat 
exchange between the furnace and the charge on the shape of the growth 
interface. 

The two-dimensional concentric fin model assumes that the axial 
location of the growth interface is the same within both the inner and 
outer charge fins. The error in predicted radial temperature variations 
obtained through the use of this assumption is expected to increase with 
interface curvature but presently is not known quantitatively. A two- 
dimensional computer model as described above may be used to determine this 
error. 

1)16 accuracy of a computer model of Bridgman growth as described above 
will be compromised by thermal prc^rty data of questionable reliability. 
Unfortunately, the thermal properties of m 2 iny semiconductors, especially 
those of the liquid, are not accurately known. In such cases, data of 
increased reliability must be found by experimental means. 

Adverse interface curvature is expected for the growth of 
semiconductor materials due to radial temperature gradients created at the 
interface. Section 8.2.1 offers several methods to counteract this effect. 
Detailed evaluation of these (and any other) methods must take place in 
order to solve this potentially serious problem inherent to the Bridgman 
growth technique. 

The technique presently considered most promising for counteracting 
the interface effect is locally heating the crucible near the interface in 
conjuncticxi with a thin crucible. Analysis of these conditions is required 
to examine the potential of this method. The relationship between the 
thickness of the crucible and the amount and distribution of heat transfer 
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to the surface of the crucible to the shape of the growth interface will 
require the solution of an appropriate two-dimensional thermal model of the 
interface region. 

This work considers furnace boundary conditions as an independent 
variable whichf in combination with the other system parameters, determines 
the interface shape. It would be useful, however, to be able to solve the 
inverse problem: i.e., what must the furnace boundary conditions be in 
order to produce a desired interface shape? There may be multiple 
solutions or no solution at all. In the first case, such results would 
provide a direct indication of the thermal design required for satisfactory 
interface shape; in the second case, the knowledge that no set of furnace 
boundary conditions exist which can produce the desired interface shape 
indicates that other parameters of the growth system require alteration. 
Ihe ability to make such conclusiois will require the computer solution of 
the appropriately modeled inverse heat transfer problem. 


- 213 - 


references 


CHiG^NAL FACE fS 

OE POOR QUALITY 


1. Brice, J. C., "The Growth of Crystals from Liquids", North-Holland 
Publishing Co., Amsterdam, London-American Elsevier Publishing Co. 
Inc., N.Y., (1973). 

2. Gilman, J. J., (ed.), "The Art and Science of Growing Crystals", Wiley, 
N.Y., (1963). 

3. Gatos, H. C., "Crystal Growth: Available Methods", Proc. Third Int. 

Summer School for Cryst. Growth, Uhiv. of New Hampshire, (1977). 

4. Flemings, M. C., "Solidification Processing", McGraw-Hill, New York, 
(1974) . 

5. Jackson, K. A., and Hunt, J. D., Acta Met. H, 1212 (1965). 

6. Schaefer, R. J., and Glicksman, M. E., Met. Trans, i, 1973 (197D) . 

7. Sharp, R. M., and Bcllawell, A., J. Crystal Growth i, 334 (1970). 

8. Tiller, w. A., Jackson, K. A., Rutter, J. w., and Chalmers, B., Acta 
ret. i, 428 (1953) . 

9. Riquet, J. P. , and Durand, F., J. Crystal Growth ii, 303 (1976). 

10. Clyne, T. W., J. Crystal Growth 691 (1980). 

11. Cline, T. W., J. Crystal Growth 684 (1980). 

12. Coriell, S. R. , and Sekerka, R. F., J. Crystal Growth ifi, 479 (1979). 

13. Coriell, S. R., Boisvert, R, F., Rehm, R. G., and Sekerka, R. F., J. 
Crystal Growth M, 167 (1981) . 

14. Chang, C. J., "A Study of Natural Convection and Its Effects on Melt 
Crystal Growth", PhJ), Thesis, Massachusetts Institute of Technology, 
Dept, of Qiem. Eng., Sept. 1982, 

15. Chang and Brown, in preparation for J. Crystal Growth. 

16. Stockbarger, D., C. , Review of Scientific Instr. 2, 133 (1936). 

17. Chang, C. E., and Wilcox, W. R., J. Crystal Growth 21, 135 (1974). 

18. Davis, K. G., Can. Met. Quart. 11, (2) 317 (1972). 

19. Bartholemew, D. M. L., and Hellawell, A., J. Crystal Growth 453 
(1980) . 

20. Vetter, J., Fruehauf, J., and Schneider, H. G., Crystal Res. and Tech, 
li, (12) 1401 (1981). 


- 214 - 



ORIGINAL PAGc IS 
OF POOR Ql’'-' *TY 


21. Fur T.-W., and Wilcox, w. R., J. Crystal Growth 51, 557 (1981). 

22. Sukanek, P. C. , J. Crystal Growth 5&/ 208 (1982). 

23. Sukanek, P. C., J. Crystal Growth 5fi, 219 (1982). 

24. Sen, S., and Wilcox, W. R., J. Crystal Growth 2&, 36 (1975). 

25. El-Mahallawy, N. A., and Farag, M. M., J. Crystal Growth Mf 251 
(1978) . 

26. Fu, T.-W., and Wilcox, W. R., J. Crystal Growth ifi, 416 (1980). 

27. Domanski, R., Fiszdon, J. F., Styczek, A., ASME paper no. 81-WA/HT-39, 
(1981) . 

28. Jones, C. L., Capper, P., and Gosney, J. J., J. Crystal Growth 5fi, 581 
(1982) . 

29. Carslaw, H. S., and Jaeger, J. C., "Conduction of Heat in Solids", 2nd 
ed., Oxford (diversity Press, Amen House, Lcxidon (1959). 

30. Rohsenow, W. M., and Choi, H. Y., "Heat, Mass and Momentum Transfer", 
Prentice-Hall, 1st ed., Englewood Cliffs, New Jersey (1961). 

31. Arpaci, V. S., "Conduction of Heat in Solids", 1st ed., Addison-Wesley, 
Reading, Massachusetts (1966) . 

32. CRC Handbook of Chemistry and Physics, The Chemical Rubber Company, 
Cleveland, Ohio, 51st ed. (1970-1971) . 

33. Eckart, E. R. G., and Drake, R. M., "Analysis of Heat and Mass 
Transfer", 1st ed., Mcgraw-Hill, New York (1972). 

34. Devienne, F. M., "Low Density Heat Transfer", in Advances in Heat 

Transfer, v.2, Acadenic Press, New York, 1965. 

35. Kraus, A. D., "Analysis and Evaluation of Extended Surface Thermal 
Systems", 1st ed.. Hemisphere Publishing Corp., Washington D. C. 
(1982) . 

36. Boyce, W. E., and DiPrima, R. C., "Elementary Differential Equations 
and Boundary Value Problans", 2nd ed., Wiley, N.Y, (1969). 

37. Ogata, K., "Modern Control Engineering", Prentice-Hall, Englewood 
Cliffs, N.J., (1970). 


- 215 - 


APPENDIX A 

DERIVATIQN OF BQUATIQN 4.11 




<Cic fS 


OF POOR QUALITY 


This ajipendix presents the derivation of eq. [4,11], 

In nondimensional form, the heat conduction equation in the charge is 
given by eq. C4.12]. Integrate each term in eq. [4.12] according to the 
formula: 



0/2 

\2— 


Q(r,2)Z.-n-y' 



Ar 





[A.1] 


(Subscript "m" in eq. (A.l] denotes an area-averaged quahtity.) 
Equation [4.12] becomes: 



[A.2] 


The boundary condition at the surface of the charge is, in nondimensional 
form: 


IM-\ 




[A.3] 


where 63 is the surface temperature of the charge. Substituting eq. [A.3] 
into eq. iA.2] yields: 


12 


a 

AC" 




AC. 


= 0 


[A.4] 


The crucible effect (section 4.2.2) is included by using the effective Biot 
and Peclet numbers in eq. [A.4] resulting in eq. [4.11]. 
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Accordi..g to the moving fin model described in Chapter 4, eq. [4,17] 
gives the axial temperature distribution in each of the four regions listed 

below: 


1. Hot zone. 

2. Liquid part of the gradient zone. 

3. Solid part of the gradient zone. 

4. Cold zone. 


The solution of eq. [4.17] for each zone is; 
where: n » 1» 2r 2, 4 


[B.l] 


(The subscript "n" denotes the zone number according to the list above). 
The constant term, 1.0, in eq. [B.1] is the particular solution for the hot 
zone and is zero for the other zones. The characteristic roots, W4 r are 

given by: 




w,n 







[B.l] 


Since the moving fin model assumes that Bic - 0, the characteristic roots 
in the gradient zone can be simplified: 
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/. 

II 

0 

' Pel 


[B.3] 

W.'s * 0 




There are eight unknown coeffici^ts in eq. two for each of the 

four zones. For the results given in this appendix, it is assumed that the 
length of charge in the hot and cold zones is infinite; under these 
conditions, B] > A 4 > 0. Boundary conditions between adjacent zones, 
given by eqs. [4.18], [4.19] and [4.20], are used to calculate the 

remaining coefficients. The boundary condition of equality of temperature 
between adjacent zones, eq. [4.18] yields: 

I2. 

“ Aj 1 5 IB.4I 

where: j/. ^ 

« axial location of the growth interface 

At the ends of the gradient zone, the boundary conditions of equality of 
heat flux between adjacent zones, eq. [4JL9] , gives; 



“ '^2 ^ 


-w 





[B.5] 
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The flux boundary condition at the interface, eq. [4.201, yields: 


. A. 


'^2 

e. * 





where: 



1 + Ks(S‘-0 


* _ Rw 


[B.6I 


Eqs. [B.4], [B.5] , and [B.61 present six linear equations in the six 
unknown coefficients A^, A 2 , B 2 , A 3 , B 3 , and B 4 . After solving for these 
coefficients, the axial gradient in the liquid near the interface, and 
the ncndiinensional interface temperature are given by: 



[B.7] 


[B.81 


(Secti(xi 4.4.3 explains that the axial gradient at the interface calculated 
from the one-dimensional model must be interpreted as the average of the 


- 219 - 


ORIGINAL PAGE IS 
OF POOR QUALITY 

charge and cctxible axial gradi^its. d&iotes this average value.) Note 
that eq. [B.81 requires an iterative solution if it is desired to determine 
5; given 

For small Pe*, eq. (5.6] holds and WjJ ^4 can be approximated 

by eq. (5.71. This leads to a simplified form of the above ^ressions: 


-z*Ice,cry'^x<.-2g,'lPe,'R; 





1 + R* I 


■"‘-Xa-ZCil 

i 

1 + 




The assumptiais that BIq ■ 0 and yield relatively simple 
analytical results. When it is desired to r?lax these assumptions, the 
preceding development can be appropriately modified. When Bi(; ft 0, for 
example, Wm^2 ^ '^m,3 longer described by eq. [B*31. Additionally, 
a particular solution must be included in eq. [B.l] for the gradient zone 
regiais; the form of the particular solution depends on the temperature 
distribution in the gradient zone annulus. When and/or Xt is not 
infinity, the coefficients B^ and/or A 4 is not zero. Additional boundary 
conditions at the ends of the charge must then be employed. In these cases 
it is more convenient to determine the coefficients and B^ by computer. 
This has been the procedure used to determine the results presented in 
section 5.5 for finite charge length. 
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CALCULATTQN OF HEAT TRANSFER QDEFFICIE^^S 

Heat transfer between the furnace walls and the outer crucible surface 
occurs by conduction and radiatiorj across the annular gap which separates 
them. The gap width is here considered sufficiwtly small so that natural 
convective flow of the gas in the gap does not occur. Ihe magnitude of the 
conduction and radiation heat transfer is here described by heat transfer 
coefficients, h^ and hj respectively. Since the conduction and radiation 
heat transfer occurs in parallel, the overall h is simply the so;, '<f h^ and 
hr = 

(If - ' k (Tf - T„) ic.li 


quality 


SAMPLE HEAT TRANSFER CDEFFICIEMS 

For an anni'lai' gap whose width is small compared to the diameter, the 
conduction heat transfer coefficient can be approximated by: 



[C.2] 


Table C.l gives heat transfer coefficients across a 1 mm gap for air and 
helium at a number of temperatures. (Helium has a high thermal 
conductivity compared to most gases.) For different gap widths, eq. [C23 
can be used to adjust the sample values given in the table. 

Radiation heat transfer across the gap, as shown in Fia. C.l, is 
described by the following equations assuming that the surfaces are grey; 

^'V(f Tocr)= IC.41 
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crucible (cr) 


C.l: GeoTietr'/ used to descri-be radiation heat transfer 

bet'-een the furnace and the outer crucible surface. 
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Teirperature 

CO 

Helium 

Air 

kgap ^32] 

he 

kgap 133] 

he 

500 

0.00299 

299 

0.000564 

56.4 

700 

0.0035 

355 

0.000663 

66.3 

900 

0.00409 

409 

0.000769 

76.9 

1100 

0.00462 

462 

0.000876 

87.6 


he measured in W/m^-K 
kgap measured in W/cm^K 
tgap - 1 nin 


Table C.l: Conduction heat transfer coefficients across a 

1 mn thick gas layer. 



v^ere; to cr) » radiation heat flux from the furnace to 

the outer crucible surface (c.f., 

Fig. C.l) per unit surface area of the 
crucible 

cr » 5.729 X 10-8 w/m2-K^ 

D » diameter 
€ ■ surface emissivity 
Tf » furnace temperature 
Tcj- = crucible temperature 
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If the gap width is small compared to the diameter of the charge (i.e., 
Df fC.4] becomes; 

T = 1 

_l_ J_ j [C.5] 

Further, if the surfaces are black, ; in this case, the 

radiation heat transfer is maximized as “ i* 

When Tf is not extremely different from Ter/ tC.3] can be 

conveniently approximated in order to yield a radiation heat transfer 
coefficient in sirapl form [301; 


where; T = (Tf + Tcr-)/2 

Equation [Cl?] shows that the radiation heat transfer coefficient increases 
as the absolute temperature cubed. Table C.2 lists sample values of the 
radiation heat transfer coefficient for various T, calculated according to 
eq. [C71 and assuming that =1. 

RAREFIED GAS QDNDUCTIQN 

The values listed in Table C.l for the conduction heat transfer 
coefficient assume that the pressure of the gas in the gap is sufficiently 
large so that the gas can be considered as a continuum. As the gas 
pressure decreases to a sufficiently low value, however, the gas no longer 
behaves as a continuum and the conduction heat transfer across the gap 
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T(C) 

hp (W/m2-K) 

500 

106 

700 

211 

900 

370 

1100 

593 
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Table C.2; Radiation heat transfer coefficients between black 
surfaces which are separated by a narrow auinular gap. 


decreases. Devienne [34] distinguishes the following limiting regimes for 
conduction heat transfer in a gas; 


< 0.001 
isfe- > 10 


continuum regime [C.8a] 


free molecule regime [C.Sbl 


where: t^p » mean free path loigth of a molecule 

of the gas 


m the free molecule regime, a gas molecule has little probability of 
colliding with other gas molecules between collisions with the furnace and 
crucible walls; in the continuum regime, this probability is very high. 

The mean free path length of a gas molecule, Ljnfp, is given by [30]; 


Lwfp = 0. ^0 ? / cr n 


[C.9] 
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where; 


<r » 

d ■ "effective" diameter of the gas molecule 
n » number of molecules per unit volume 


The ideal gas law may be used to express the molecule density of the 
gas, n: 


r, > 1C.101 

^ R^T 

where; N » Avagadro's number, 6.02 x lO^^/gm-mole 
» universal gas constant 
T » temperature 
p ■ pressure 

Subsituting eq. [C.lOl into [C.9] the mean free path length becomes; 


L 



<Tp NJ 


[C.ll] 


Helium and air are again used as typical examples. Their effective 
diameters, d, are [30]; 


He; d » 2.18 x 10“8 cm [C.12a] 

Air; d * 3.72 x 10“® cm [C.12b] 

Using these values, eq. [C.ll] gives the mean free path length at a 
temperature of 900 *C as; 


He; 



N - CCW^) 


[C.13a] 
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p Ctorr) 


Cfm 


EC. 13b] 


Equations EC.131 can be used in eqs. tC.8] in order to determine the 
pressure limits for the continuum and free molecule regimes. 

In the free molecule regime, the conduction heat transfer coefficient 
is given by Deviaine [34]: 



where: y = Cj/Cy 

Cp « specific heat at constant pressure 
Cy ■ specific heat at constant volume 
M = molecular weight of the gas 
a « accomodation coefficient 


EC.14] 


The parameter a in eq. EC.14] depends on the effectiveness of energy 
transfer during collisions between a gas molecule and the walls of the 
furnace and crucible. If this energy transfer is complete, the value of a 
is equal to unity and heat transfer is maximized. A value of unity is 
presently used so that eq. EC.14] will give an upper limit for the free 
molecule ccxiduction heat transfer coefficient. A reasonable approximation 
to the temperature T in eq. EC.14] is the mean value of the furnace and 
crucible surface temperatures. Evaluating eq. EC.14] at T*900 *C, for 
helium and air: 

He: h^EW/m^-K) = 47.6 p(torr) EC. 15a] 

at 900 “C 

Air: hcEW/m^-K) = 34.6 p(torr) EC.lSb] 
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At pressures required so that conduction heat transfer across the gap is in 
the free molecule regime, eq. [C.15] shows that the conduction heat 
transfer coefficient is much smaller than the radiation heat transfer 
coefficients given in Table C.2. Conduction heat transfer is, in this 
case, essentially eliminated. 


TEMPERATURE 




5 OF THE HEAT TRANSFER (D EFFICIEOT 


The thermal models developed in Chapters 4 and 6 assume that the Biot 
number is not a function of position within the hot, gradient and cold 
zones. The temperature dependence of the conductive and radiative heat 
transfer coefficients, as indicated in Tables C.1 and C2, affects the use 
of these models as follows: 


(1) The temperature distribution in the system is not known 
a priori. Since the heat transfer coefficients are 
temperature dependent, the Biot numbers are also not 
initially known. An estimate of the temperature levels 
must be made in order to calculate an initial estimate of 
the heat transfer coefficients. From the resulting 
temperature response of the system the accuracy of the 
initial estimate can be checked. New estimates can be 
made for the heat transfer coefficients and the procedure 
is repeated until the desired level of accuracy is 
obtained. 

(2) Since the temperatures are not spatially uniform within 
the system, the heat transfer coefficients are also not 
spatially uniform. A reasonable estimate of a single 
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appendix q 


END BOUNDARY QQNDTTTONS AND INFINITE FURNACE LENGTH 


nj. END BOUNDARY CDNDITIQN 

When the length of the charge in the hot and/or cold zone is not 
infinite (see section 5,5), boundary conditions at the end of the charge 
must be included in the thermal model of the charge. The hot end of the 
charge, for example, receives heat from the hot zone furnace by radiation 
and convection (Fig. D.l). A heat transfer coefficient describing the heat 
transfer between the hot end of the charge and the hot zone furnacb can be 
approximated in order to calculate a Biot number for use in the hot end 
boundary condition. If Ah > , the magnitude of this Biot number has 
little effect on conditions near the interface; on the other hand, if 
Ah h ' value of the hot eid Biot number is required in order 
to completely specify the thermal model of the charge at the hot end. 

The cold end of the charge is attached to the lowering rod and, hence, 
the cold end boundary condition is more complicated. A simple thermal 
model for this case, shown in Fig. D.2, is presently used in order to 
calculate an approximate Biot number for the cold end boundary condition. 
A pull rod of thermal ccxiductivity kp, diameter Dp, cross sectional area Ap 
and length Lp, is attached to the cold end of the charge; the heat transfer 
between the furnace and the pull rod is described by a heat transfer 
coefficient, hp. The pull rod and charge are assumed to be in perfect 
thermal 'contact, the cold zone heat pipe is assumed to be infinitely long, 
and radial temperature gradients are neglected. Motion of the system is 
neglected (i.e., Pe is assumec negligible). 
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The cold end of the charge experiences a heat loss to the pull rod 
which may be used in order to define a heat transfer coefficient for the 
cold end . ^ the charge, h^,p: 

- - kf Ay ».n 

where: " axial distance from the charge/pull cod boundary, 

^ see Fig. D.2 

(The subscript "cp" denotes the charge/pull rod boundary.) 
Nondimensionalizing eq. [D.l] defines a cold aid Biot number for the cold 
end boundary condition: 

Sic, le.. (3^ - O) - ,D.2a, 

g' B ^ I [D.2bl 

’’ ■ k» A iCf, 

The thermal conditions at the ends of the pull rod are described by the 
"transmission" matrix developed by Kraus [3b] . In nondimens ional form, the 
transmission matrix equation for the pull rod is; 


M m 




■> m 

~q) 


1 

1 

tm “ 






L acp J 


[D.3] 
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■ nondiinensional axial coordinate, based on diameter Dp, 



■ effective Biot number for the pull rod, 

■ Biot number for the pull rod used to calculate Bip ; 


The transmission matrix relates the temperature and the gradient in 
temperature at one end of the pull rod to these conditions at the other 
end. Equation tD.31 may be used to develop an expression for the axial 
gradient of temperature in the pull rod at the charge/pull rod boundary: 


If the pull rod is very long, the second term within the brackets of 
eq, [D.4] becomes negligible; eq. [D.4] then takes the form: 






Substituting eq. [D.5] into eq. [D.21 : 


[D.5] 



[D.61 


Bi given by eq. [D.61 can be used as an approximation for the cold end 
cp 

boundary condition of the charge, i.e., eq- [D.2bl. 

Equation [D.6] gives the Biot number for the cold end boundary 
condition ior an infinite length pull rod of uniform properties and cross 
section. The transmission matrix concept 135] can also be conveniently 
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employed to approximate the end Biot number for pull reds of varying 
composition and cross section. Equation [D.6] also applies for extensions 
at the hot end of the charge. 

The charge would appear infinitely long at on of its ends if it was 

attached to an extension at that end which has the same thermal properties^ 

cross sectional area and heat transfer coefficient as the charge itself. 

In this cassf eq. [D.6] indicates that the end Biot number is twice the 

square root of the local charge Biot number: 2 Jb!^. This leads to the 

conclusion, reported by Chang and Wilcox [17], that the charge appears 

I — 

infinitely long if the end Biot number is 2jBij^^^. Bence, eq. [D.6] 
indicates the manner which the end Biot number may be adjusted, by the 
proper choice of kp, Ap and/or Bi**, in order to approximate infinite 
length conditions with an extenstion attached to the end of the charge: 

INFINITE FURNACE LE?K5TH 

An assumption in the thermal model of the Bridgman, system io that the 
furnace temperature boundary condition is of infinite extent in both the 
hot and cold ends (section 2.2.2). The actual case of finite furnace 
length is illustrated in Fig, D.3 for the cold zone furnace. The part of 
the charge not within the furnace experiences different boundary conditicwis 
than that part of the charge within the furnace. The following brief 
analysis explains under what conditions the furnace boundary conditions may 
be considered infinite in extent, thereby making unnecessary the inclusion 
of the thermal conditicxis outside of the furnace in the thermal models. 
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The analogous expression to eq. [D.4] for the charge within the cold 
zone of Fig. [D.3] is: 

[D.8] 

where: 

Equation [D.8] indicates that the thermal ccxiditions at the hot end of the 
cold zone (i.e.f at ) are little affected by the conditions at 

the end of the cold zone furnace if: 

acXCw^^^Xc.') 1 

Equation [D.9] is analogous to the infinite charge length criteria 
developed in section 5.5. Therefore, if the length of the furnace 
satisfies eq. [5.14], the furnace lo:agth can be considered infinite. This 
conclusion is valid for either the hot or cold zone furnace. 
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TBANSiani: 




AL BEHAVIOR 


The dimensionless form of eq. C4.1] for the one-dimensional thermal 
model of the charge is: 


^ - ?e g"* - 9 „) - = 0 

where: X ® 

The nondimensicnal time, t , in eq. CE.1] is measured with reference to the 
time required for the charge to be lowered a distance equal to one charge 
diameter . 

Uie axial coordinate in the present work, C , is measured with respect 
to a stationary reference frame (i.e., the center of the gradient zone); 
the transient term of eq. [E.l] therefore accounts for changes in the 
charge temperature with time at a given location within the furnace (i.e., 
not the changes in temperature of a material portion of the charge). The 
quasi-steady assumption in the thermal model of section 2.2.2 implies that 
the transieit term of eq. [E.1] is negligible. 

^ THERMAL TRANSIENTS CMJSED BY NDN-INFINITE CHARGE LENGffl 

Temperature transients do not occur when the charge length is infinite; 
therefore, when the length of charge within the hot and cold zones 
satisfies the infinite length criteria, eq. [5.14], the quasi-steady state 
assumption is valid. 

When the charge length can not be considered infinite, temperature 
transients in the charge occur. Results obtained from a quasi-steady model 
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may still, howav-^r, provid# inf ormatj or-, which is useful in a qualitative 
sense. For example, the curves of lie. 5.8 indicate that the axial 
position of the interface moves from the cold to the hot furnace as the 
charge is lowered, Riquet and Durand 1 91 and Sukanek 1 22, 231 also neglect 
transient heat transfer in their studies of the dependence of the growth 
rate on the axial position of the char e 

Large Biot numbers decrease thi infinite charge length criteria, 
eg. [5.141; a larger portion ot the charge is solidified under quasi-steady 
conditions as the Biot number increases. An alternative to larger Biot 
numbers for making the charge uppear infinitely long is to increase the 
apparent length of the charge t x)th the hot and cold ends by the addition 
of extensions to the crucible. vSee Appendix D.) This is normally 
accomplished at the cold end of the charge through its thermal 
communication with the pull rod. 
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Section E.l discusses temperature transients caused by the steady 
motion of a finite length charge tlirough the furnace, when a sudden change 
in one of the furnace operational j-\arameters ve.g., lowering rate, hot and 
cold :one temperatures) occurs, the coefficients of eg. 1E.1] (i.e., Pe and 
likewise change .and temperature transients occur. The neglect ot the 
transient term under these circumstances would indicate a suixien shift from 
the guasi-steady temperature field in the charge before the change to the 
new quasi-stead%* state after tlie chance. 

Fu and Wilcox [211 liave studied the transient interface pi\sitiv.vi after 
a sudden change in the lower invi rate for a oharoe of intinite lenotli. 
Their results indicate that the interface chances with time from its 
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initial to final quasi-steady position in an exponential fashion; 





tE.2] 


where: i “ quasi-steady interface position after change in 

■’ lowering rate 

« quasi-steady interface position before change in 
lowering rate 


The time constant of the exponential function in eq. [E.2 ] , ^ , is shown to 
decrease for larger Biot numbers, smaller gradient zone length, Lq, and 
smaller latent heat factor, Rr. They suggest that the time constant is 
similar for transients caused by a sudden change in furnace temperature. 
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D BEAPMEWr QF ms FIN SQIIATIQNS FOR mE TO D -DIMENSIONAL THERMAL MODEL 

With the assumptions of Sections 2.2 and 6.2, each of the 
concentrically located fins (i.e., inner charge, outer charge, crucible, 
gradient zone annulus) can be considered as a moving thin fin. According 
to Carslaw and Jaegar [29J , the axial temperature distribution for the fins 
in the hot and cold zones is given by the following equations: 








0 


\ocAout ^ 2 ^ 


"Tout Tc,r 


oot Hootj 


= 0 


cr 






LJSL - a ^ ^ 

.0 


[F.la] 


[F.lb] 


[F.lc] 


where; A » cross sectional area of a fin 

R » thermal resistance per unit length between fins 
subscripts; 

loc » local charge phase, liquid or solid 
in ■ inner charge fin 
out ■ outer charge fin 
cr ■ crucible 

f » hot or cold zone furnace 


When considering the fins within the gradient zone, Tf and Rcr,f in 
eq. [F.lc] change to Tq and Rcr,G, where the subscript "G" denotes the 
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The radial geometry of the concentric fins shown in Fig. $.2 
specifies the cross section areas and thermal resistances in eqs. tF.l]. 
The cross section areas Ain are chosen equal; therefore, 
Din » D/JJl The radial locations for the representative temperatures of 
the fins are Di/2, D 2/2 and D 3/2 for the inner and outer charge fins emd 
crucible fin respectively; they are chosen so that there is an equal area 
within the fin to either side; 


D, = t)/2L 

[F.2a] 

D^ = JV4 T> 

[F.2b] 


[F.2c] 


The thermal resistances can then be expressed as (Rohsenow and Qioi [26]): 


K 


IrtjOW I 


2-^ fei,c 


[F.3a] 



[F.3b] 








- s 


-I 




2 tt K 


or 




[F.3c] 
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h ■ heat transfer coefficient from furnace or gradient 
zone annulus to the surface of the crucible based on 
the crucible outer surface area 


Substituting eqs. [F,2] into eqs. into eqs. tP.l] and non-dimensionalizing 
results in eqs. [6.1] and [6.2]. 
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SION FOR APPROXIMATE 


RFACE Q3RVRTURE 


This appendix employs simple approximations in order to develop a 
useful relation between isotherm curvature and the axial and radial 
temperature gradients obtained from the one- and two-dimensional thermal 
models. 

An isotherm within the charge is assumed to be spherical; its 
curvature can then be represented by a single number — its radius. The 
cucial distance between the locations of an isotherm at the center aind the 
surface of the charge is denoted by The assumed isotherm geometry, 
shown in Fig. G.l, yields the following relation between and the 
isotherm curvature: 


T) Z 4N ^ 

where N is the radius of curvature of the isotherm in number of charge 
radii. (The reciprocal of N is the curvature.) Eq. [G.l] gives a positive 
value of N when the isotherm is curved as shown in Fig. G.l. If the axial 
gradient does not vary greatly over the cross section: 


^ ^sc 

where: '4) - 4’if^o) 

Substituting eq. [G.21 into [G.l] ; 


[G.2] 


[G.3] 
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For the purpose of relating <p^ , the heat conduction equation 

within the charge is used: 




f ''' 



[G.41 


If the Pe term is neglected and if it is assumed that is not a 
function of ■j , the resulting temperature distribution in the radial 
direction is parabolic. In this case# it is easily shown that: 



Using eq. [G.5] in eq. [G.3] : 


[G.5] 




[G.6] 


Equation [G.6] indicates that isotherm curvature is small (i.e., large N) 
when the axial temperature gradient is large as well as when the radial 
temperature gradients are small. 
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